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On  the  Design  of  a  Towed  ELF  H-Field  Antenna 
^  I.  Introduction 

Over  the  last  few  years,  the  possibility  of  building  an  ELF  loop 
antenna  (H-field  antenna)  for  towing  behind  a  submarine  has  been  receiving 
close  attention.  The  principal  result  of  this  work  is  the  conclusion  that  the 
task  can  be  accomplished  orifplWien  the  many  problems  are  thoroughly 
understood  and  then  accounted  for  in  the  antenna  design.  Although  at  pre- 
,>-  *8 ent  no  design  of  sufficiently  good  performance  exists,  much  of  the  basic 
understanding  has  been  obtained.  Thus  the  work  from  now  on  should  con¬ 
sist  largely  of  using  the  understanding  as  a  basis  for  testing  and  evaluating 
candidate  construction  materials  and  fabrication  methods.  To  facilitate 
the  work,  this  report  documents  that  understanding.  It  does  not,  however, 
propose  «  specific  antenna  design  or  predict  eventual  performance.  This 
cannot  be  done  until  a'ter  the  testing  and  evaluating  stage. 

Interest  centers  on  the  frequency  range  of  about  20  to  200  Hz  in  the 
ELF  band  because  tn  that  range  the  attenuation  in  the  water  and  in  the  earth 
.exosphere  waveguide  is  low  enough  for  long-range  broadcasting  to  sub¬ 
merged  submarines  to  be  ieasible  [6],  The  problem  of  radiating  enough 
energy  from  an  anteruia  of  practical  size  is  a  difficult  one,  but  it  appears  to 
be  solvable  by  r  ing  a  horizontal  electric  dipole.  This  kind  of  antenna  can 
also  bo  u.ec  as  the  receiving  antenna  on  the  submarine.  Physically  it  con- 
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s ists  of  two  small  electrodes  built  into  a  cable  and  separated  by  some  200 
to  300  m  and  towed  well  aft  away  from  the  interfering  electromagnetic  noise 
fi  'id  generated  by  the  boat  (Fi&.  1).  It  confers  on  the  boat  an  ability  to 
communicate  over  a  very  useful  range  of  depths  and  speeds.  Unfortunately, 
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it  is  linearly  polarized  and  the  downward  propagating  transmitted  signal  in 
the  ocean  is  also  mainly  or  wholly  "linearly  polarized.  It  suffers,  therefore, 
from  the  disadvantage  that  it  becomes  .nearly  or  wholly  blind  if  the  submar¬ 
ine  heading  is  perpendicular  to  the  direction  vector  of  the  transmitted 
electric  field,  in  the  far-zone  field  of  the  transmitting  antenna  this  occurs 
when  the  submarine  heading  is  perpendicular  to  the  great-circle  path  linking 
the  transmitter  and  receiver. 

The  successful  development  of  an  H-field  antenna  built  into  the 
same  cable  would  provide  a  receiving  polarization  perpendicular  to  that  of 
the  E-field  (eiec.trode-pair)  antenna.  Since  the  transmitted  electric  and 
magnetic  field?  are  mutually  perpendicular  and  essentially  horizontal,  the 
H-t’ield  antenna,  in  the  form  of  a  long  slender  solenoidal  vinding,  would  be 
maximally  sensitive  when  the  E-field  antenna  is  blind.  Communication  can 
therefore  be  maintained  independently  of  heading  (Fig,  2). 

ihe  fundamental  limitation  on  ELF  communication  at  the  surface  is 
atmospheric  noise.  Compared  with  it,  the  noise  generated  by  the  E-field 
antenna  and  the  interference  generated  by  the  boat  are  both  negligibly  small. 


a* 


Fig.  1.  Electrode -pair  E-field  antenna  towed  from  the  submarine. 


Fig.  2.  Combined  E-  and  H-field  antennas  for  polarization  diversity. 
(Ce  and  are  factors  dependent  upon  the  geometrical  details  of  the 
two  antennas.  ) 


As  the  antenna  submerges,  both  signal  and  atmospheric  noise  are  attenuated 
by  the  same  amount  and  so  the  signal  to  noise  ratio  is  initially  not  affected. 
The  communication  rate  is  therefore  unimpaired.  However,  the  local 
noises --from  the  antenna  and  the  boat- -are  essentially  unaltered  by  sub¬ 
merging.  If  the  antenna  depth  continues  tc  increase,  eventually  the  atmos- 
pheric  noise  level  will  be  equal  to  that  of  the  local  noises,  and  communica¬ 
tion  begins  to  be  degraded.  Thus  the  depth  of  operation  is  determined  by 
the  antenna  noise.  The  design  problem,  therefore,  for  the  H-field  antenna 
is  to  reduce  the  level  of  the  noise  it  generates  to  be  no  greater  than  that  of 
the  E- field  antenna. 

It  is  shown  in  Section  III  that  the  thermal  noise  generated  by  an 
H-field  antenna  without  a  ferromagnetic  core,  given  the  limitations  of 
available  conducting  materials  and  the  constraint  that  the  antenna  cable  be 
positively  buoyant,  is  too  high  unless  the  antenna  is  allowed  to  become  in- 
conveniently  large.  One  is  forced,  therefore,  to  consider  using  a  ferro¬ 
magnetic  core.  Unfortunately,  t'ns  gives  rise  to  "magnetostrictive" 

(stress -induced)  noise  and  Barkhausen  noise.  The  former  arises  because 
the  geomagnetic  bias  flux  density  in  the  core  changes  its  magnitude  when 
the  core  is  stressed  by  mechanical  vibration,  thereby  inducing  a  change  in 
voltage  in  the  antenna  winding.  The  latter  occurs  because  the  geomagnetic 
bias  flux  density  in  the  core  does  not  follow  smoothly  the  changing  axial 
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component  of  the  geomagnetic  field  strength  as  the  boat  turns.  Rather,  it 
proceeds  by  a  successioif  of  small  jumps,  which  induce  voltage  transients 
in  the  antenna  winding,  [l,  pp.  613  and  524], 


Magnetostrictive  noise  can  be  reduced  by  choosing  the  right  core 
material,  by  keeping  the  net  bias  flux  density  close  to  zero  by  means  of  an 
adjustable  direct  current  in  the  antenna  winding  [4],  by  making  the  antenna 
long  and  less  sensitive  towards  the  two  ends  than  it  is  in  the  middle  of  its 
length,  by  minimizing  the  stress  transmitted  to  the  core  and  by  reducing  the" 
level  of  the  cable  vibration. 

Barkhausen  noise  can  be  reduced  by  choosing  the  right  core  mat¬ 
erial,  by  keeping  the  net  bias  flux  density  as  constant  as  possible  [4]  and 
by  increasing  the  length  of  the  antenna. 

Noise  is  also  generated  when  the  geomagnetic  flux  linking  the  turns 


of  the  antenna  winding  is  modulated  by  their  angular  vibration.  This  is 


called  motion- induced  noise  and  is  present  whether  or  not  a  ferromagnetic 
core  is  used.  Reducing  the  cable  vibration,  making  the  antenna  long, 
making  it  leas  sensitive  towards  the  two  ends  that  it  is  in  the  middle  of  its 
length  and  making  its  variation  of  sensitivity  along  its  length  very  smooth 
are  techniques  for  reducing  this  noise.  A  comparison  of  the  motion-induced 
noise  of  an  H-field  antenna,  when  this  is  due  to  ideal  cable  flexing,  with 


that  of  an  E-f.eld  antenna  shows  that  in  principle  a  long  H-field  antenna  need 
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be  no  more  susceptible  to  motion-induced  noise  than  the  E-iield  antenna  [7]. 


Recently,  however,  it  has  been  point eji  out  that  more  complicated  cable 


motion  is  possible  than  that  due  to  ideal  flexing,  but  no  quantitative  esti¬ 


mate  is  made  of  tne  relative  magnitude  of  the  resulting  noise  [193. 


there  is  as  yet  no  evidence  that  such  motion  is  an  important  noise  source. 


On  the  other  hand,  the  point  is  valid  in  principle  and  deserves  closer 


attention. 


T.  date,  two  Tong  H-tield  antennas  have  been  made  and  tested. 


Neither  is  quiet  enough  to  be  useful,  but  their  performance  taken  together 


with  the  improvements  likely  to  be  realized  by  applying  the  additional  tech 


niques  covered  in  this  report,  suggest  that  a  usefully  quiet  antenna  can  be 


built.  One  antenna  showed  that  a  low  enough  thermal  noise  can  be 


attained,  although  it  suffered  from  uncontrolled  motio'n- induced  noise  due 


to  independent  vibrations  of  its  core  elements,  while  the  other  antenna 


showed  that  a  core  construction  exists --the  helically  wound  tape  core-- 


which  eliminates  the  offending  type  of  core  vibration,  although  it  suffered 


from  a  high  thermal  noise  and  a  vib-ation  induced  noise  that  appears  to  be 


ascribable  to  its  uniform  sensitivity  di'tribution  (see  Section  VII), 


!>mn 


The  following  Sections  treat  in  detail  the  points  mentioned  in  this 
brief  summary,  except  for  the  newly  purposed  motion-induced  noise  mechan¬ 
ism  mentioned  above  which  has  yet  to  be  fully  explored.  A  brief  survey  of 
the  experimental  evidence  is  also  included.  It  gives  strong  support  to  those 
aspec  theory  Vv  touches. 
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II,  Measure  of  Performance  and  Effective  Length 

A  convenient  measure  of  antenna  quality  is  one  that  allows 
antennas  of  different  configurations  and  types  to  be  compared  and  which, 
at  the  same  time,  uses  the  minimum  number  of  additional  assumptions.  It 
should  also  be  some  physically  identifiable  quantity,  if  possible,  to  aid  in 
its  application. 

A  measure  that  fits  this  prescription  well  is  the  ratio  of  the 
r.  m.  s.  open-circuit  noise  voltage  to  the  effective  length  of  the  antenna, 
where  the  effective  length  is  the  length  by  which  an  incident  electric  field 
must  be  multiplied  to  give  the  corresponding  open-circuit  antenna  voltage. 
The  measure  has  the  dimensions  of  volts  per  meter  and  is  physically  the 
horizontal  electric  field,  evaluated  at  the  depth  of  the  antenna,  which  pro¬ 
duces  an  antenna  voltage  equal  to  the  noise.  In  general,  the  ocean  conduc¬ 
tivity  and,  for  antennas  with  an  appreciable  vertical  dimension,  a  depth 
reference  point  on  the  antenna  must  be  specified.  No  assumptions  about 
the  strength  of  the  transmitted  field  or  of  the  atmospheric  noise  field  are 
required.  Neither  is  it  necessary  to  specify  any  particular  depth  of 
operation. 

The  term  "equivalent  noise  field"  or  ENF  will  be  used  to  denote 
this  measure  of  performance. 

For  the  trailing  E-field  (electrode-pair)  antenna,  the  effective 
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length  is  simply  the  distance  between  the  electrodes.  Thus  the  total  ENF 
for  this  antenna  is  its  measured  noise  voltage„d\vided  by  the  electrode 
separation.  Typical  ENF  spectra  at  8  and  14  knots  are  shown  in  Fig.  3 
for  a  300-meter  electrode-pair  antenna  T3].  For  this  length  of  antenna, 
the  noise  is  now  believed  to  be  dominated  by  motion-induced  noise. 

Under  the  conditions  14  knots  and  45  Hz  (to*be, adopted  here  as  the 
standard),  Fig.  3  shows  the  ENF  of  the  300  meter  E-field  antenna  to  be 
-205  dBE. 

For  a  horizontal  H- field  antenna  immersed  in  the  ocean,  the 
effective  length  can  be  calculated  using  Faraday's  law.  If  H  is  the  axial 
component  of  the  incident  magnetic  field,  then  the  induced  voltage  per 
meter  length  of  antenna  is  given  by  u>|icHAN,  where  u>  is  the  radian  fre¬ 
quency,  |i  the  reversible  permeability  of  the^ynVenna  core,  A  the  cross- 
sectional  area  of  the  core  and  N  the  turns  density  (turns  per  meter)  in  the 
antenna  winding.  The  total  antenna  voltage  V  is  given  therefore  by 

V  =  '»hJ'\ic  ANdz  =  f  e  E  (1) 

where  the  integral  is  carried  out  over  the  complete  antenna  length,  and  E 
is  the  incident  electric  field  associated  with  the  magnetic  field. 

Both  signal  and  atmospheric  noise  propagate  down  into  the  ocean 
from  the  surface.  At  ELF,  the  field  components  E  and  H  are  essentially 
horizontal  and  mutually  orthogonal.  They  are  related  by  the  wave  imped¬ 
ance  of  the  ocean,  V  -iu>p0/ <7  ,  o  being  the  conductivity  of  the  water  and 
the  permeability  of  free  space.  Using  this  relation  in  (1),  one  obtains, 


u  zu  60  80  100  120  MO  160  180 

FREQUENCY  (Hi) 

Fig.  3.  Typical  equivalent  noise  field  spectra  for  a  300-meter  trailed 
E-field  antenna  at  two  towing  speeds.  (From  Griffiths  [3]. ) 


(2) 


in  magnitude 

ie  =  T  /  KaANdz 

where  pa  =  pc/pQ  *s  the  relative  small-signal  permeability  of  the  core  and 
6,  defined  as  V  2/ (too  pQ)  ,  is  the  *kin  depth  for  electromagnetic  propagation 
into  the  ocean. 

*  V 

This  formula  has  been  obtained  by  neglecting  the  magnetic  flux 
in  the  area  occupied  by  the  antenna  winding.  This  flux  density  links  the 
outer  layers  of  the  antenna  winding.  The  permeability  of  this  part  of  the 
cross  section  is  that  of  free  space,  however.  Thus  (2)  is  valid  provided  the 
addition  of  the  core  increases  the  antenna  inductance  by  a  large  factor. 

The  effective  length  of  an  H-field  antenna  can  be  accurately 
estimated  by  measuring  its  impedance  in  aix.  or  in  the  ocean.  Appendix  A 
shows  that  its  impedance  in  the  ocean  is  negligibly  different  from  its  im¬ 
pedance  in  air,  and  the  inductance  Lq  in  air  is  given  simply  by 


AN  dz  • 


(3) 


Thus,  if  the  turns  density  is  uniform,  (2)  and  (3)  can  be  combined  to  give 


_  V?L{ 


6  (i  N 
o 


(4) 


If  the  turns  density  is  not  uniform,  but  the  pYodtic t  pcA  is,  then  (2)  and  (3) 


yield 


f  = 
e 


1 


N  dz 


lA, 


N2  dz 


(5) 


in  which  the  ratio  of  integrals  can  be  evaluated  from  the  winding  specifica¬ 
tion  of  the  antenna. 
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III.  Thermal  Noise 


The  thermal  noise  of  the  antenna  is  the  noise  generated  by  the 
resistance  of  the  antenna  measured  at  the  antenna  terminals.  This  resist¬ 
ance  is  the  sum  of  the  resistance  of  the  connecting  wires,  the  resistance 
of  the  antenna  winding  and  the  "radiation"  resistance  of  the  winding,  which 
represents  the  losses  in  the  antenna  core  and  in  the  surrounding  ocean. 

The  core  losses  can  be  calculated  using  standard  methods  [2,  p.  78  ff.  ] 
and  can,  in  any  case,  be  avoided  by  subdividing  the  core  to  interrupt  the 
induced  circulating  currents.  The  losses  in  the  water  are  difficult  to  cal¬ 
culate  in  general  but  Appendix  A  shows  that  their  contribution  to  the  ENF 
is  very  small.  The  resistance  of  the  connecting  wires  can  be  made  small 
compared  with  that  of  the  winding  by  using  a  large  enough  number  of  turns 
on  the  antenna  winding.  Thus  the  critical  component  of  the  thermal  noise 
is  that  contributed  by  the  antenna  winding. 

It  is  instructive  to  examine  'what  can  be  achieved  without 

a  ferromagnetic  core  in  the  antenna.  In  Appendix  B,  it  is  shown  that  even 
if  all  the  weight  of  the  cable  were  devoted  to  the  winding  and  if  insulation, 
where  necessary,  could  be  infinitesimally  thin,  then  the  minimum  thermal 
ENF  attainable  in  a  300  m  length  of  the  current  0.65  in.  O.  D.  buoyant 
antenna  cable  at  45  Ha  is  about  -180dBE.  This  could  be  reduced  to  the 
required  -205  dBE  by  increasing  the  diameter  by  a  factor  of  about  4  or 
increasing  the  length  by  a  factor  of  about  300.  The  latter  is  certainly  im¬ 
practicable;  the  former,  since  it  implies  an  increase  of  volume  by  a  fac¬ 
tor  of  sixteen,  is  still  excessive.  A  ferromagnetic  core,  th<  refore,  would 
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seem  to  be  essential. 


Since  the  reversible  permeabilities  attained  by  currently 
available  ferromagnetic  materials  can  be  several  tens  of  thousands  times 
larger  than  that  of  free  space,  the  effective  length  of  the  antenna  can  be 
increased  enormously,  without  any  increase  in  the  thermal  noise.  Thus 
the  thermal  ENF  can  be  reduced  substantially.  The  difficulties  are  that 
i)  a  long,  slender,  high  permeability  core  may  be  magnetically  saturated 
by  the  geomagnetic  field  and  thereby  lose,  temporarily,  its  high  reversible 
permeability,  ii)  the  stretching  and  flexing  of  the  cable  may  plastically 
strain  the  core  material,  causing  a  permanent  loss  of  its  high  reversible 
permeability,  iii)  the  longitudinal  cable  vibrations,  transmitted  to  the  core, 
generate  magnetostrictive  (i.  e.  ,  strain-induced)  noise,  and  iv)  the  varia¬ 
tions  of  the  axial  geomagnetic  field  component  due  to  heading  changes  can 
generate  Barkhausen  noise. 

Geomagnetic  core  saturation  can  be  prevented  by  using  a  core 
material  of  relatively  low  permeability  or  by  driving  a  direct  current  of  the 
appropriate  magnitude  and  sense  through  the  antenna  winding  to  bring  the 
net  flux  density  in  the  core  close  to  zero.  This  latter  technique  can  be 
made  automatic  by  using  an  out-of-band  sensing  signal  to  monitor  continu¬ 
ously  the  inductance  of  the  antenna  and  make  appropriate  changes  in  the 
bias  current  [4]. 

The  plastic  strain  problem  must  be  dealt  with  by  careful 
mechanical  design,  and  techniques  for  reducing  magnetostrictive  noise 
and  Barkhausen  noise  are  discussed  in  later  sections. 

The  addition  of  a  te rromagnetic  core  raises  the  question  of 
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determining  the  optimum  proportions  of  winding  and  core.  The  buoyancy- 

restriction  places  a  limit  on  the  total  (reduced)  weight  per  unit  length  of 

the  core  plus  winding.  Thus  the  minimum  thermal  ENF  is  achieved  at  a 

proportion  somewhere  between  the  all-conductor  one.  with  its  low  effective 

length  and  the  all-core  one,  with  its  very  high  resistance. 

If  a  is  the  core  radius,  b  the  outer  winding  radius,  2,0 the  antenna 

length  and  a_  the  wire  conductivity, then  assuming  that  (j.  ,  the  average  rela- 
w  a 

tive  permeability  of  the  core,  is  large  compared  with  unity,  the  effective 
length  of  the  antenna  is  4 

a  -  ZvaZ,x  >jz  h£./6  (6) 

G  d 

while  the  total  winding  resistance  is 

„  _  2irN20  b+a 
CTw  b  -  a 

assuming  in  each  case  I’-  it  the  turns  density  N,  the  permeability,  pa>  and  the 
geometry  are  independent  of  length.  Thus,  from  (6)  and  (7),  the  power 
spectral  density  Se (to )  of  the  thermal  ENF  is  given  by 


S  (to) 
e 


kb  Tk  s‘ 


b  +  a 


2  4 

£  a  (b"a)- 


(8) 


where  is  Boltzmann's  constant  and  T^  the  absolute  temperature  (°K). 

On  the  other  hand,  the  reduced  mass  per  unit  length  w*  of  the 
antenna  is  given  by 


i  2  1  ,,  2  2 ,  | 

w 1  =  ir  a  P  +  ir  (b  -  a  )  0,  ' 

C  w 


(91 


where  P  '  is  P  -0.  and  p  '  is  P  -  P.  .  The  densities  P  ,  0  .  and  P,  are 
c  (  i  w  w  b  c  w  b 

those  of  the  tore,  the  winding  and  the  buoyant  flotation  material  (foamed 
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polyethylene  in  the  current  buoyant  antenna  cable).  The  reduced  densities 
Pj  and  pj  appear  in  (9)  because  the  antenna  core  and  winding  displace  not 
free  space  or  air  but  the  buoyant  flotation  material. 

The  problem,  then,  is  to  minimize  (8)  while  keeping  (9)  fixed 
at  the  maximum  allowed  by  the  buoyancy  restriction.  This  is  easily  done, 
■Tor  at  this  optimum  state,  the  total  differential  of  (8)  and  of  (9)  must  both 
be  zero.  That  is 


dS  (u>)  =  0  = 


2kbV 
17  Va2* 


- -  C(2a2+ab-2b2)da  -a2dbl 

a^b-ar 


dw'  2  0  =  2Tr[(Pc' -  P^)ada  +  Pjbdb]  . 

For  these  equations  to  have  a  solution,  the  determinant  of  their  coefficient 

matrix  must  he  zero.  That  is 

(2a2  +  ab  -  ab2)  P  'b  a3(p  '  -  P  ')  =  0 
W  c  W 

or,  with  3  =  pc'/pw  =  ^cTV^w''  Pb*  and  a  =  b/a* 

2<y3  -  <y2  -  2<y  -  8  +  1  =  0.  (10) 

The  functional  relationship  between  a  and  8  is  presented  graphically  in 
Fig.  4. 

Thus,  given  the  densities  Pc>  Pw  and  P^,  one  calculates  8  and 
then  solves  (10)  (or  uses  Fig.  4)  to  evaluate  the  optimum  radius  ratio  b/a. 
From  this,  the  individual  radii  a  and  b  can  be  found  by  using  (9)  and 
these  then  substituted  in  (8)  to  give  the  minimum  value  of  the  thermal  ENF 
power  spectrum. 
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It  is  not  necessary  in  practice  to  achieve  precisely  the  optimum 
proportions,  because  large  deviations  from  the  optimum  cause  relatively 
minor  degradation  in  the  thermal  ENF.  To  demonstrate  this,  one  notes 

that  ifw'  and  O  '  are  given,  then  (9)  shows  that  l/a  is  proportional  to 

^0*  2  -**' 

0  +<*  -  1,  which  means,  from  (8),  that  the  mean  square  thermal  ENF  is 

2  2 

proporiional  to  (0  i  or  -  1)  (a  +  l)/(o-  1).  This  function  is  plotted,  with  a 

or 

(equal  to  b/a)  as  the  abscissa,,  for  various  values  of  0  in  Fig.  5.  The 
^.curves  show  that  the  relative  winding  thickness  a-  1  ,<sa^  vary  over  a  range 
of  at  least  three  to  one  without  increasing  the  minimum  thermal  ENF  by 
more  than  1  dB. 

The  total  reduced  mass  per  unit  length  that  the  cable  can  sup- 

port  is  given  simply  by  tr  c  (P  -  p^),  where  c  is  the  outside  cable  radi,us  and 

D  in  the  average  density  required  for  the  complete  cable.  However,  some 

of  this  (say,  w&')  may  have  to  be  allocated  to  strength  members  or  some 

other  essential  piece  of  the  structure.  Thus  the  available  reduced  mass 

2 

per  length  w1  to  be  allocated  to  the  antenna  is  ttc  (p  -  p^)  -wa'. 

In  the  current  buoya/il  antenna  cable,  for  example,  c  is  0.  325 
inches,  is  0.63pQ  and  the  18  fiberglass  strength  members  are  each 
0.  038  inches  in.  diameter  and  have  a  density  of  about  1.  67  P  .  (Here  Pq 
is  10  kg/m  ,  the  mass  density  of  water).  Thus  w.'  is  18  x  (0.  019x  0.  0254) 

a 

x  it  x  (1. 67  -  0.63)  and  so,  if  P  is  to  be  0.  SSp^,  then  w'  =  0.  0334  kg/m. 

But,  if  the  core  density  Pc  is  8.  5  pq  and  the  winding  density  Pw 

is  2.7,  then  0  =  3.802  and,  from  (10),  b/a  =  1.636.  Using  (9),  one  then 

-  3 

finds  that  a  =  0.  968  x  10  m.  taking  6  to  be  37,  5m,  k  to  be 


1.  38  x  10"23  j/°K,  Tk  to  be  300°K,  <Jwto  bo  3.57X  107  mho/Wand  21  to  be 
300  m,  one  finds  that  iW'drde.r  fop  the  the mml  ENF  to  be  -205  dBE,  the 
relative  core  permeability  must  be  746, 

it  • 

In  practice,  the  crosB-floctlonnl  of  the  ferromagnetic  part 

of  the  core  will  be  some  fraction  y  of  the  total  core  croBB  soction.  If  the 

remaining  space  is  filled  with  a  material  having  tho  density  PR,  then  the 

average  density  of  this  composite  core  Ib  given  by  p  e  Vp*  +  (l  -  Y)p  and 

c  c  n 

its  average  relative  permeability  la  given  by  |i&  =  whejt4o'pc  and  are 
the  density  and  relative  p^/neability  of  the  ferromagnetic  component  of  the 
composite  core.  Table  I  shows  the  result  of  calculating  the  relative  perme¬ 
ability  required  for  the  thermal  ENF  to  be  -205  dBE  for  a  range  of  values 
of  the  core  utilization  factor  Y  »  The  assumptions  are  that  p*c  =  8.  5Pq  and 
p„  =  o  .  The  core  diameter  is  a  function  of  Y  because  P  depends  on  Y  and 
hence  the  optimum  proportions  vary  with  Y. 


Required  u  fo*  *he  Thermal  ENF  to  be  -205  dBE 


Y 

Core  O.  D. 
(inch) 

b/  a 

f*r 

.  005 

0.291 

1.  0814 

24794 

.  01 

0.  280 

1.  0877 

13000 

.  02 

0.260 

1.  0998 

7073 

.05 

0.  220 

1.  1333 

3450 

.  1 

0. 182 

1.  1820 

2174 

.  2 

0.  144 

1.  2622 

1471 

.  5 

0.  101 

1.4364 

966 

1.  0 

0.076 

1.  6355 

746 

Table  I 


T 


w  ■ 


■  • 


The  relative  permeabilities  specified  in  Table  I  apply  to  an 
antenna  having  a  uniform  sensitivity  profile.  A  more  likely  profile  is  one 
that  is  parabolic,  because  it  discriminates  better  against  vibration  noise 
(both  magnetostrictive  and  motion-induced).  In  this  case,  the  permeability 
would  have  to  be  raised  by  the  factor  3/2  to  maintain  the  same  antenna  ef¬ 
fective  length.  And,  if  the  ferromagnetic  material  is  applied  in  ihe  form 
of  a  helix  of  helix  angle  9,  then  the  permeability  must  be  increased  by  the 
further  factor  l/cos^9. 

It  should  be  noted  that  the  numbers  appearing  in  Table  I  are  the 

result  of  a  particular  set  of  assumptions.  Changing”the  assumptions  can 

cause  a  large  change  in  the  required  permeability.  For  example,  if  P  (the 
average  mass  density  of  the  .complete  antenna  cable)  is  taken  to  be  pQ  instead 

of  0.  85  pQ,  thenw’  is  nearly  twice  as  large,  which  means  that  the  required 

permeabilities  are  only  one  quarter  of  the  values  shown. 

The  permeabilities  shown  in  Table  I  are  mostly  small  compared 
with  the  high  permeabilities  attained  by  modern  ferromagnetic  alloys.  The 
mechanical  constraints  imposed  by  this  application  are  severe,  however. 
They  make  it  impossible  to  realize  even  a  small  fraction  of  this  potential. 

As  a  last  resort,  the  cable  diameter  could  be  increased,  for,  as  (8)  shows, 
if  all  transverse  dimens-’ons  are  doubled,  the  required  permeability  is  re¬ 
duced  by  a  factor  of  four. 

IV.  Motion  -  Induced  Noise 

If  v  (z,t)  is  the  transverse  displacement  of  the  antenna  center 
line,  as  a  function  of  time  and  of  position  along  the  antenna,  from  its 
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nominally  straight  condition,  then  by/bz  is  the  tangent  of  the  angle  it  makes 
with  the  z  axis.  Provided  |by/bz|<j<  1,  therefore,  the  geomagnetic  flux 
linkage  per  unit  length  is  given  by  A(z)N(z)(Bz  +  By  Sy/bz),  where  A(zV  is  the 
area  of  the  antenna  winding,  assumed  for  the  moment  to  have  an  air  core, 
and  N(z)  is  its  turns  density  (turns  per  meter).  The  geomagnetic  field  is 
assumed  to  lie  in  the  yz  plane,  so  only  its  two  components  By  and  Bz  appear. 

The  total  moticn  -induced  voltage  v(t)  is  given  by  the  integral 
over  the  length  of  the  antenna  of  the  time  rate  of  change  of  the  flux  linkage 


per  unit  length.  Thus,  since  3Bz/dt  =  3By/bt  =  0, 
v(t)  =  B^AWNI.)  ez. 


(HI 


In  Appendix  C  it  is  shown  that  the  same  equation  applies  if  the 
antenna  has  a  ferromagnetic  core,  except  that  A(z)N(z)  must  include  the 
additional  factor  p  ,  the  relative  incremental  permeability  averaged  over 
the  core  area. 

This  expression  can  be  simplified  in  the  case  of  the  uniform 
solenoidal  winding,  for  then  A(z'Nfz)  is  a  constant,  say  AqNo,  over  the 
length  of  the  antenna,  axA  is  zero  elsewhere.  The  expression  then  reduces 
to 

v(t)  =  ByAoNQ  —  {y(!,t)  -  y  (-!,«}  ,  (12) 

b  t 

assuming  that  the  antenna  extends  from  -£  to  +£  along  the  z  axis.  Thus  only 

the  vibration  of  the  two  end-points  generates  the  noise  voltage,  no  matter 

1 

how  long  the  antenna  is.  By  simply  preventing  the  ends  from  vibrating, 
one  can  reduce  the  noise  voltage  to  zero  even  though  the  rest  of  the  antenna 
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vibrates  in  a  completely  arbitrary  way. 

It  is  impractical  to  prevent  a  particul'a-r^oint  on  the  towed  cable 
from  vibrating,  but  this  uniform -winding  example  immediately  suggests 
another  way  of  reducing  the  motion-induced  noise  voltage.  For  the  noise 
voltage  is  clearly  associated  with  'e  discontinuities  at  ±i  of  the  "sensitivity 
..  profile"  A(z)N(z).  Where  the  profile  is  uniform,  no  noise  is  generated. 

Thus  a  profile  which  is  as  uniform  as  possible  (within  the  constraint  that  it 
must  be  zero  for  jz|>f)  should  generate  the  least  noise. 

One  cannot  simply  remove  the  discontfcSBiiies  at  the  ends  of  a 
uniform  profile.  For  then,  being  zero  at  ±f,  the  profile  must  be  zero 
everywhere  and  the  antenna  no  longer  exists.  One  seeks  therefore  a  non- 
uniform  profile  which  goes  to  zero  at  ±1.  For  this  class  of  function,  an 
integration  by  parts  converts  (11)  to  the  form 


where  the  more  compact  notation  u(z)  has  been  used  for  the  sensitivity 
profile  A(z)N(z)  and  u'(z)  stands  for  du(z)/dz. 

It  appears  that  in  the  frequency  band  of  interest  here,  the  sources 
of  the  transverse  cable  vibrations  are  the  stress  fluctuations  at  the  cable 
surface  due  to  the  turbulent  fluid  flow  past  the  obstacle  presented  by  the  cable 
itself  (see  Appendix  D),  and  the  large  attenuation  for  transverse  mechanical 
waves  along  the  cable  (see  Appendix  I)  isolates  the  antenna  from  disturb¬ 
ances  generated  bv  the  ship.  Thus  the  motion-induced  noise  arises  from  a 
statistically  uniform  distribution  of  transverse  forces  which  are  essentially 


uncorrelated  from  point  to  point  along  the  caMe.  And  since  they  are  uncor¬ 


related,  the  total  voltage  is  equal  to  the  square  root  of  the  sum  of  the 
squares  of  the  voltages  generated  by  each  elemental  force  separately. 

To  calculate  the  separate  voltages  it  is  more  convenient  to  work 
in  the  frequency  domain,  in  which  (13)  can  be  .expressed  as 

V(u>)  =  «.By/ u'(»)Y(z,  o>)dz  ,  (14) 

where  V  and  Y  are  the  Fourier  transforms  of  v  and  y.  But  in  this  notation, 
the  displacement  Y  (z,iu)  caused  by  a  single  transverse  force  F(io)  applied 
at  the  point  zq  is  given  by  G(z,  zq,  tu)  F(u>),  where  G(z,  zq,  u>)  is  the  Green's 
function,  and  is  the  displacement  at  z  due  to  a  unit  force  of  radian  frequency 
u)  applied  at  zq.  Hence  the  voltage  induced  by  this  force  is,  from  (14), 

V(oj)  =  iu)  By  j'u'  (z)G(z,  zq,  u>)dz  F(oj)  ,  (15) 

The  power  spectral  density  ^v(u>)  of  the  antenna  voltage  induced  by  this 
single  force  is,  therefore,  from  (15),  given  by 

§v(u>)  =  u)2  By^  Sj. (u>)  )  /u'(z)C(z,  V„)dzl  (16) 

where  Sj.(oj)  is  the  power  spectral  density  of  the  force.  Thus,  if  there  is 
a  uniform  distribution  of  M  of  these  forces  per  unit  length,  each  one  s'atis- 
tically  independent  of  the  others,  then  the  total  power  spectral  density 
S  (u>)  of  the  antenna  voltage  is 

Sv(u>)  =  -i)2  By2  MSf  (uj)  J I /  u'(z)  G(z,  zo,  cu)dz  |  dzQ  .  (17) 

Now  MS^((u)  is  related  to  S  j.(k,  u>),  the  two-dimensional  power 
spectral  density  of  the  transverse  force  distribution  on  the  cable,  defined  as 
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sf(k,  to)  =  ycf(c,  T)e“ikC+iU,TdCdT  (18) 

where  C,(C,t),  the  two-dimensional  correlation  function,  is  the  expected 
value  of  f(z,t)f(z+£,  t  +  t )  and  f(z,  t)  is  the  transverse  force  per  unit  length. 
Since  the  mean  square  total  force  on  a  piece  of  cable  which  is  'long  com¬ 
pared  with  the  correlation  length  of  the  force  distribution  is  8^.(0, tu)  times 

A 

the  length,  on  the  one, band,  and  MSf(u>)  times  the  length,  on  the  other,  then,~^' 

A 

the  relationship  is  simply  MSj(id)  =  Sj.(0,  u>). 

The  antenna  voltage  power  spectrum  can  therefore  be  rewritten 


as  C  C  2 

Sv(u>)  =  co?'  By2  Sf (0,  co )  l\  I  u'(z)  G(z,  zQ,U))dz  |  dzQ  ,  (19) 

Further  analysis  is  complicated  by  the  fact  that  the  tension  in 
a  towed  cable  increases  linearly  from  the  free  end  towards  the  tow  point. 
This,  in  turn,  means  that  the  mechanical  dynamics  of  the  cable  are  non- 
uniform,  and  so,  in  particular,  that  G(z,  zq,co)  depends  upon  the  sum  of 
the  pair  z,  z  in  addition  to  their  difference.  One  fundamental  result  can 
be  obtained  fairly  simply,  however,  on  the  assumption  that  the  sensitivity 
profile  u(z)  is  both  perfectly  parabolic  and  also  very  long  compared 
to  cne  linear  extent  of  the  region  of  influence  of  a  force  applied  at  a  point 
ri.e.  ,  long  compared  with  the  "width"  of  G(z,  z  ,  co)].  If  these  conditions 
are  met,  the  u'(z)  is  essentially  a  constant  over  the  "width"  of  G(z,  zq,  to ) 
and  so  can  be  taken  outside  the  inner  integral  and  assigned  the  value  u'fz  ). 


But  the 


n  /  G(z,  z 

J  ° 


,  co)dz,  assuming  that  the  tension  is  essentially  invariant 


over  the  "width"  of  G(z,  zQ,  co),  can  be  calculated  by  assuming  the  tension 
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is  constant  at  the  value  it  takes  at  zq.  Appendix  E  shows  that  the  integral 
is  then  given  by  - l/(to  m  ),  where  m.  is  the  effective  transverse  vibrating 
mass  per  unit  length  of  cable.  Thus,  for  u(z)  sufficiently  long  and  suffic¬ 
iently  smooth,  (19)  can.be  rewritten  as 


B 


Sv(o»  = 


:27-V-sf<° 


®  lmt'‘ 


,  .)  ^ I  u,(zq)  I  dzQ 


(20) 


The  motion-induced  ENF  power  sp'ectrumv#5'  (to)  is  given,  therefore,  by 
dividing  (20)  by  the  effective  length  squared.  Therefore,  from  (2)  and  (20), 
and  using  the  relation  u(z)  =  u  (z)  A(z)  N(z),  one  finds 

cL 


For  a  given  cable  size,  cable  buoyancy,  measuring  frequency 
and  towing  speed,  all  quantities  in  the  right  member  of  (21)  are  predeter¬ 
mined  except  for  u(z).  Tlius  the  problem  is  to  choose  u(z)  to  minimize  the 
auotient  of  integrals  appearing  in  (21)  with  the  constraint  that  u(z)  is  con¬ 
tinuous  for  all  z  and  zero  for  lz!  By  applying  the  calculus  of  variations 

one  finds  the  solution  to  be 


uo(z) 


uo(0)[i-(z/f)2],  |z!  <  e 

o  ,  U.\*« 


(22) 


and  this  will  be  referred  to  as  the  "parabolic  sensitivity  profite"  on  account 
of  its  shape,  shown  in  Fig.  6.  For  this  optimum  profile,  the  quotient  of 
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Fig.  6.  Optimum  sensitivity  profile  (parabolic  profile)  for  minimmng 
the  motion -induced  ENF. 


3 

integrals  reduces  to  the  expression  3/(2f  ),  and  so 


S  (u>)  =  -=rj- •  f  y  )  St(0,  0)) 
6  4(i  l»|mt|  I  f' 


Significant  features  of  this  formula  are,  first,  that  the  tension, 
bending  stiffness  and  energy  absorption  of  the  cable  do  not  appear  in  it,  only 
the  transverse  force  spectrum  and  the  effective  mass.  Thus  since  the  buoy¬ 
ancy  is  essentially  predetermined,  the  motion-induced  ENF  at  a  particular 
speed  and  frequency  depends  only  upon  cable  diameter  and  antenna  length 
(provided  the  "long,  smooth  profile"  assumption  is  not  violated). 

Also  of  interest  is  the  fact  that  the  absolute  antenna  sensitivity 
does  not  appear  in  the  formula.  This  means  that  the  motion-induced  ENF 
depends  only  upon  the  shape  of  the  sensitivity  profile  and  not  upon  whether 
or  not  a  ferromagnetic  core  is  used,  for  example. 


In  Appendix  G,  it  is  shown  that(tl7c?5i%?/ion-induced  ENF  of  an 
E-field  antenna  (the  electrode-pair  antenna)  is  given  by  B^S^O,  uj  )/(2fu/'  |  m^  I  ^). 
Comparing  this  expression  with  (23),  one  sees  that  if  the  two  antennas  are  of 
equal  length,  then  they  generate  the  same  motion  -  induced  ENF  when  they 
are  V6  6  long.  When  they  are  shorter  than  this,  the  H-field  antenna  is 
noisier  than  the  ’E-field  antenna,  and  when  they  are  longer  than  this,  'fre 
reverse  is  true.  In  a  4  mho/m  ocean  at  45  Hz,  6  is  37.  5  m.  Thus  V*6  6 
is  91.8  m  or  about  300  feet.  (This  comparisorTtfepends  upon  the  additiorai 
reasonable  assumption  that  S^O,  <u)  is  the  same  for  both  the  x-tand  the  y- 
directed  transverse  force  distribution.  That  is,  that  the  transverse 
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buffeting  forces  are  statistically  azimuthally  uniform). 

A  more  difficult  question  to  treat  is  the  modification  of  (23) 
required  to  take  into  account  the  effect  of  a  short,  rough  sensitiv¬ 
ity  profile.  This  will  be  carried  out  under  the  assumption  that  the  tension 
is  constant  over  the  whole  length  of  the  antenna.  The  further  slight  modi¬ 
fication  required  to  account  for  varyitvg^^sion  is  examined  in  Appendix  H. 

The  basic  motion-induced  noise  equation  is  (19)  in  which,  since 
now  the  tension  is  being  assumed  constant,  G(z,  zq,  <o)  can  be  written*cfs 

G  (z-z  ,  oj).  Then,  with  the  definitions 
o  o 


U(k) 
u(z)  = 


Ht(k.u>) 


(z)  e~*kz  (jjj 


U(k)eikzdk 


Go(z>U5)  =  Zir 


~  j  H(k,U))ei 


(24) 


(u)  e"ikz  dr. 

X 


dk  , 


the  integral  J’u.'iz)  Gq(z  -  zq,  tu)  dz  can  be  rewritten  as  i(2ir)  y*kU(k)Ht(k,  to) 


exp  (ikz  )  dk  and  so  (19)  reduces  to 
o 


/2  2 

k2  |  U(k)l  |Ht(k,  <0)|  dk  .  (25) 


If  tl^e  assumption  that  the  transverse  buffeting  forces  are  essentially  un- 
correlated  from  point-to-point  along  the  cable  is  not  valid,  then  it 
follows  that  Sf(0,  m)  in  {co)  must  be  included  under  integral  sign  as 
Sj.(k,  iu).  Thus  the  "uncorrelated  forces"  assumption  is  equivalent  to  making 
the  assumption  that  the  spatial  spectrum  of  the  force  distribution  is  white 
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(uniform)  over  the  range  of  k  for  which  the  product  |kU(k)Ht(k,  to)  |  has  non- 
negligible  magnitude.  In  Appendix  D,  this  assumption  is  shown  to  be 
satisfied  in  practice. 

The  functions  |U(k)|  and  |H (k,u>)J  appearing  in  the  integrand 

.**&**’  '  2 

in  (25)  are  sketched  in  Fig.  7.  The  peak  in  jU(k)l  reflects  the  fact  that  to 

receive  the  signal,  the  winding  must  be  unidirectional,  so  that  u(z)  has  a 

large  "d.  c.  "  component  at  k=0.  The  peaks  of  |  Ht(k,  tu )  I  occur  at  k  =  ±  kfc  , 

whererfA‘t  is  the  real  part  of  the  ordinary  characteristic  wave  number  for 

transverse  wave  propagation  at  radian  frequency  u>  along  the  cable.  (Since 

the  cable  is  not  perfectly  flexible,  there  also  exists  an  "extraordinary" 

characteristic  wave  number.  It  has  a  large  imaginary  part,  however,  and 

2 

so  does  not  have  a  readily  visible  effect  on  the  shape  of  |  Ht(k,  o> )  1  along 
the  real  k  axis.  ) 

If  the  sensitivity  profile  u(z)  is  very  long  and  very  smooth,  then 
its  Fourier  transform  U(k)  will  be  very  peaked  at  the  origin  and  be  very 
small  elsewhere.  Then  the  integrand  in  (25)  (the  product  with  k  of  the 
two  functions  sketched  in  Fig.  7),  will  be  of  negligible  magnitude  for  all  k 
except  for  two  narrow  regions  symmetrically  placed  either  side  of- and  close 

i  t2 

to  the  point  k=0.  But  over  these  regions,  |Ht(k,u>)|  is  essentially  equal  to 

,  ,2 

|  Hfc(0 , to )  1  and  so  the  integral  in  (25)  can  be  rewritten  as 


|  Ht(0,ui)l  2/l  kU(k)|  dk  ,  which,  by  virtue  of  Parseval's  theorem,  can  be 

expressed  as  2tt^  Hfc(0 , )  I  u'(z)]  dz  .  Substituting  this  back  in  (25),  and 

2  1 

noting,  from  Appendix  E,  that  H^O,^)  ~  -l/(a>  m  ),  one  finds  that  (25)  re- 

i  i2 

duces  to  (20).  Thus  the  contribution  of  the  peak  in  |U(k)'  at  the  origin  to 


re- 
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Fig.  7.  Sketch 
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the  integral  in  (25)  is  simply  the  fundamental  long,  smooth  profile  contribu¬ 
tion  discussed  earlier. 

There  are  two  other  peaks  in  the  integrand  in  (25),  however, 

whose  contribution  must  also  be  included/'  If  the  profit.e^yfc)  is  short 

enough,  but  still  smooth,  then  the  minor  peaks  of  U(k)  will  make  the 

contribution  of  the  integrand  arising  from  the  peaks  of  H^(k,o) 

at  k  =  ±  kt  larger  than  that  arising  from  the  peak  in  U(k)  at  k=0.  Since  the 

peaks  of  H.(k,u>)  are  in  practice  very  high  and  sharp,  their  contribution  to 

2  2  f  2 

the  integral  in  (25)  can  be  written  k  Av{IU(k  )|  }  /  |H,(k,(Ju)|  dk,  where 

2  2  J 
Av{  |U(k.)l  }  is  a  local  average  of  j  U(k)|  centered  on  the  points  ±  k^.  But, 

from  Appendix  E ,J *  |Ht(k,u))!2dk  is  the  same  as  2tt  Sy(iu)/,Sj(0,  u>),  where 

Sy(o>)  is  the  spectral  density  at  frequency  to  of  the  transverse  displacement 

of  the  antenna  at  a  point. 

Thus  (25)  can  be  expressed  in  the  approximate  form 

,2 


S  (w)  =  uj2B  2S.(0,u>) 
v  y  f 


Au'(z)]  dz  2  2  S  (») 

- -  +k.  Av{|U(k.))  } 

I  4  2  t  t  a 

L  »  |xntl 


Sf(0,u>) 


(.16) 


n 

The  magnitude  of  Avi  I  j  depends  upon  two  things.  On  the 

one  hand,  if  the  profile  is  very  smooth,  then  U(k)  is  given  asymptotically 
for  large  k  by  (Lighthill,  [5]  ) 


U(k)  -  -  —  £  *u'(z  ie“ikZn  , 
k2  n  n 


(27) 


7$ 


where  Au'fz^)  is  the  jump  in  u'(zn)  at  the  point  of  discontinuity  z^  and  the 
summation  includes  all  points  of  discontinuity.  Thus  Av(  |U(k)|2}  is 


30 


approximately* 


AvC|U(k)j2}  ~  —  £  |Au«(zn)i 


On  the  other  hand,  if  the  profile  is*^-1  rough"  due,  for  example,  to  variations 
in  the  turns  density  or  in  the  me  ability,  then  |U(k)J,  for  k  not  close 

to  zero,  can  be  much  bigger  than  (27)  predicts.  For,  expressing  the  actual 


rough  profile  ur(z)  as 

u  (z)  =  u{z)[l  +  e  (z)l  >  (29) 

r 

where  e(z)  is  the  "fractional  perturbation  '  of  the  nominally  smooth  profile 

2 

u(z),  then  one  can  show  that  the  expected  value  of  |ll(k)|  is  augmented 
by  the  term  S  (k)  f  u2(z)dz.  Here  S  (k)  is  the  spectral  density  of  e(z)  and  is 
defined  as  the  Fourier  transfopjj^oi  the  correlation  function  of  e(z),  which 
in  turn  is  the  expected  value  of  e(z)e(z+£)  .  It  is  assumed  that  e(z)  has 

zero  mean  and  is  statistically  uniform  and  stationary  in  z. 

The  final  general  expression,  therefore,  for  the  motion-induced 
ENF  of  an  antenna  for  which  u(z)  is  continuous  but  u'(z)  has  discontinuities 


at  z  (n=l,  2,  .  . .  )  is  S  (u) )/ £  or 

n  v  e 


2  2  2 
u>  6  B  S 


S  >)  =  ~pr - 1 - - 

2/[u(z)dz 


.(o, uj )  r  fw 

Y7'  *Bt  a)4 


2  2 
(z)3  dz  |  Au'jz^) j  S^cu) 

I  ^  kfc2  Sj.(0,co) 


?  -  S  (0,u» 

H  k~S  (k.)  f uz(z)dz  -1 - 

‘  6  *  J  Sf  (0, u>) 


.  (30) 


If  the  profile  is  -ominally  the  ideal  one  def.ned  by  (22),  then 
J~  u(z)dz  is  4fu^(0)/3,  J'  u  •'.)]2dz  is  8  uo2(0)/(3f),  ^  j  Au'fz^)]  is 


2  2  2  2 

8  uq  (0 )/j2  and  J  u  (z)dzi3  I6uq  (0)f/l5,  sc  that  (30)  becomes 


s  (w)  =  -j 

e  4  r 


,  ?  2  ->  1  t“4m2S  (to) ' 

J_+i_k  Vs  (kfc)  - *- 

k2*  5  fc  *  fc)  Sf(0,co)  < 


By  comparing  this  with  the  simple  "long,  smooth  profile"  formula  given  by 

(23),  one  sees  that  as  the  antenna  becomes  shorter,  a  correction  term  pro- 
-4 

portional  to  f  will  eventuaVt\^be  significant,  and  that  the  effect  of  profile 
roughness  is  to  add  a  correction  term  proportional  to  the  spectral  density 
of  the  relative  roughness  and  to  f""*.  Since  the  long,  smooth  profile  term 
is  0{f  },  and  the  "shortness"  correction  is  Off  this  third  roughness 

term  will  eventually  dominate  when  the  antenna  is  made  arbitrarily  long. 

A  further  significant  difference  between  the  correction  terms 
for  shortness  and  for  roughness  on  the  one  hand,  and  the  long,  smooth 
profile  term  on  the  other,  is  that  (-he  former  are  sensitive  to  the  stiffness 


and  absorption  properties  of  the  cable.  To  take  simple  account  of  this  fact, 
it  is  convenient  to  define  a  mechanical  Q-factor  for  transverse  cable  vi¬ 
bration  as  follows 


c  1  f 

1  2  k t  J 


On  the  assumption  that  S^(k,u))  is  essentially  constant  over  the  range  of  k 

4  2 

for  which  H,(k,iu)  has  signifU~-  '  agnitude,  Q.  is  equal  to  (tr/kju)  I  m ,1  S  (to)/ 
1  'I t  c  1  t  y 

Sj.(0,<ju),  so  that  (31)  can  be  rewritten 

.  /  6B  \2  (  3  Q  2k  3f2Q. 

S  (u))  =  J-  f - S  (O.oj)  j  1  +  - L  +  - * - ~  S  (kf)  , 

413  V*KI  /  \ 


W  , 
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and  this  formula  applies  to  the  nominally  parabolic  profile.  The  is  a 
function  of  frequency  and  tension.  A  graph  of  Q^,  calculated  for  a  speci¬ 
men  of  the  standard  0.  65"  D  foamed -polyethylene  buoyant  antenna  cable, 

I*"’' 

is  given  in  Appendix  I,  together  with  a  similar  graph  of  k^. 

•**  Equation  (33)  shows  that  if  S^O.tu)  is  known,  then  by  calculating 
k  and  from  measured  mechanical  properties  of  the  cable,  one  can 
calculate  the  motion-induced  ENF  of  an  H-field  antenna  of  given  length  and 
pro/ile  roughness.  can  be  estimated  from  boundary  layer  data  (see 

Appendix  D)  or  calculated  once  and  for  all  from  measurements  of  the  vibra¬ 
tion  of  a  towed  cable  in  a  towing  tank.  Thus  (33)  frees  one  from  the 
necessity  of  carrying  out  an  expensive  and  time-consuming  submarine  test 
for  each  modification  of  the  antenna  design. 

For  the  envisioned  operational  antenna  length,  k^iis  very  much 
larger  than  (see  Appendix  I),  and  so  the  "profile-shortness"  correction 
term  in  (33)  is  small  compared  with  unity.  The  "profile- roughness" 
correction  term  can  be  very  much  greater  than  unity,  however.  Since  both 
this  Ji^Y  and  the  motion-induced  ENF  of  the  E-field  antenna  are  propor¬ 
tional  to  J l  *  (see  Appendix  G),  a  criterion  for  the  profile  smoothness  for 
the  H-field  antenna  can  be  obtained  which  is  independent  of  l.  Thus,  if  the 
prof’ le- roughness "  part  of  the  motion- induced  ENF  of  the  H-field  antenna 
is  to  be  no  greater  than  the  motion-induced  ENF  of  the  E-field  antenna,  then 
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(34) 
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This  criterion  may  be  difficult  to  satisfy  in  practice,  since  both  kfc6  and  Qj_ 
are  numerically  large. 

Since  the  bending  stiffness  and  mass  distribution  of  an  actual  cable 
are  not  absolutely  constant  but  vary  slightly  from  point-to-point  along  the 
cable,  and  since  the  "profile  roughness"  is  a  serious  noise  source,  the  ques¬ 
tion  arises  of  whether  there  should  be  a  "cable-properties-variation"  terrfi 
in  addition  to  the  two  ot’ff&Jf&or recting  terms  in  (33)-,.  It  is  easy  to  see, 
however,  that  such  a  term  is  negligible  compared  with  the  "profile  roughness" 
term.  For  if  the  profile  is  perfectly  smooth,  then  the  variation  of  the  cable 
properties  has  no  effect,  because  the  cancellation  effect  is  a  purely  geo¬ 
metrical  one  depending  only  on  the  assumption  that  the  disturbance  attenuates 
as  it  propagates  ..way  from  the  disturbing  force.  Thus  the  "cable-properties - 
variation"  term  arises  o*>V?,  h  the  profile  is  rough  and  so  represents  a  per¬ 
turbation  of  a  correction.  It  is,  therefore,  negligible. 

If  the  long,  smootv  /vofile  condition  can  be  achieved,  then  the 
advantage  of  using  a  parabolic  rather  than  a  uniform  profile  can  be  very  great, 
for,  from  (12)  and  the  definitions  of  and  (Eqs,  (2)  and  (32)),  one  deduces 

that  the  motion-induced  ENF  of  an  H-field  antenna  with  a  unifot-m  profile  is 

.2 

V(0,o>)  .  (35) 

?  \  m  I  r>1  !  /  1 

4  it 


■st>1  - 


k,°t  /  5Bv  V 
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ThuSi  from  (35)  and  (23),  one  sees  that -the  uniform  profile  generates 

noise  by  the  factor  Qfck  i/ (3tt)  than  the  purely  parabolic  profile.  For  en- 

visioned  operational  antenna  lengths  incorporated  in  the  current  design  of 

4 

antenna  cable,  this  factor  can  be  of  the  order  of  10  or ^O^dB,  as  reference 
to  Appendix  I  will  verify. 

A s  for  the  effect  of  cable  diameter  on  motion-induced  ENF, 
equation  (33),  together  with  the  expression  for  S-(0,(w)  given  in  Appendix  D, 
shows  that  the  long-smooth-profile  approximation  is  proportional  to  the 
Inverse  square  power  of  the  cable  radius  (assuming  the  cable  density  does 
not  vary).  On  the  other  hand,  the  profile-roughness  contribution  is  pro¬ 
portional  to  the  inverse  seven-halves  power  of  the  cable*radius,  if  the 
bending  stiffness  is  the  dominant  restoring  force.  For  low  frequencies  and 
large  tensions,  however,  when  the  cable  tension  provides  a  considerable 
part  of  the  restoring  force,  the  precise  law  is  difficult  to  state  simply. 
Thi^va  because  an  increase  in  radius  increases  the  bending  stiffness  and 

therefore  also  increases  the  mechanical  damping.  Thus  the  Q  is  a  compli- 

* 

cated  function  of  diameter.  In  the  extreme  case  of  very  high  tension,  since 
the  tension  itself  is  proportional  to  cable  radius,  the  profile -roughness 
contribution  to  the  motion-induced  ENF  is  proportional  to  the  inverse  half¬ 
power  of  the  cable  radius.  But,  in  practice,  the  tension  would  be  too  small, 
and  the  frequency  and  diameter  too  large,  for  this  behavior  to  be  observed. 

V.  ?A&^'netostrictive  Noise 

In  Section  I  the  manner  in  which  magnetostrictive  noise  is 
generated  svas  described.  Tae  various  techniques  for  reducing  it  (selection 
of  core  material,  using  a  degaussing  current,  stress-relieving  mechanical 
construction  and  profile  shaping)  were  also  briefly  mentioned.  In  this 
Section  the  two  latter  technioues  are  examined  in  greater  detail. 
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When  the  antenna  core  is  axially  strained,  there  is  a  concomitant 
"magnetic  strain"  in  the  axial  geomagnetically-induced  bias  flux  density  B  : 


AB  Af 

B  '  6  i 

o 


(36) 


Particular  ferromagnetic  alloys  exist  for  which  the  "strain- ratio"  e  is 

4 

minimized,  but  in  general  it  is  a  large  number  --  about  10  ,  for  example, 
for  45  Permalloy  [1,  p.  618]. 

If  the  axial  displacement  of  a  point  on  the  cable  from  its  rest 
position  be  denoted  by  x(z,  t),  then  the  axial  strain  is  6x(z,  t)/5  z  and  so  the 
change  of  flux  linkage  is  given,  from  (36),  by 

d  $  =  Bqc  NA  6x(z,  t)/Sz  , 

where  N  is  the  turns  density  of  the  antenna  winding  and  A  is  the  area  of  the 
ferromagnetic  core.  Thus  the  total  magnetostrictive  noise  voltage  is  given 
by 

/*  n 

(37) 


„n,  -  / 


e  B  N  A  dz  . 

°  &zat 


In  the  frequency  domain,  this  is 
V(tu)  =  -ill) 


)  /  e  B 

J  ° 


N A  -LX-fe»d-  dz  . 
3z 


(38) 


Now  the  axial  strain  is  present  in  the  form  of  longitudinal  me- 
chr.iical  waves  travelling  along  the  cable.  The  attenuation  of  these  waves 
u>  relatively  small  (see  Appendix  I)  and  so  disturbances  generated  outside 
the  antenna  region  travel  into  and  through  it.  Thus  X(z,u>)  has’  the  form 
X (0 , u> )  exp(ikQz)  where  kQ  is  the  (complex)  wavenumber  for  longitudinal 
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mechanical  waves.  Thus  (38)  can  be  written 

ik  z 


/US.  Zi 

e  B  N  A  e  °  dz  . 
o 


(39) 


In  practice,  the  quantities  e,  Bq  and  |JLa  are  likely  to  be  constant  along  the 
antenna.  Then,  since  u(z)  has  been  defined  in  Section  IV  to  be  the  product 


u  NAf  (39)  can  be  written  a-s 
a  ^ 


e  B 


V(u>)  =  (ok  X(0,u))  - - 

o'  \in 


U(-ko)  , 


(40) 


where  U(k)  is  the  Fourier  transform  of  u(z).  The  magnetostrictive  ENF 
spectral  density  is  therefore  given  by  dividing  V(w)  by  l  from  (2)  and 


squaring.  That  is 


se(»)  = 


■itt) 


.  *  P 


sx(.) 


ko0^ 


U(o) 


(41) 


where  S  (m)  is  the  power  spectral  density  of  x(0,  t). 

Equation  (41)  shows  that  the  antenna  sensitivity  profile  u(z) 
should  be  shaped  to  minimize  the  ratio  |U(-kQ)/U(0)  |.  Since  the  antenna 
is  likely  to  be  long  compared  with  the  wavelength  of  the  longitudinal  waves 
(i.  e.  ,  |kQf|  >:>  t*18  value  of  U(-kQ)  ig  largely  determined  by  the  magni¬ 

tude  of  the  lowest-order  discontinuity  of  the  function  u(z),  while  U(0)  is  the 
area  under  the  curve.  Mvxc5?7he  same  considerations  apply  here,  therefore, 
as  did  to  the  reduction  of  motion-induced  noise.  The  profile  should  be  is 
Jong  and  as  smooth  as  possible.  Here,  however,  the  reduction  in  noise 
obtained  by  going  from  a  uniform  to  a  tapered  profile  is  not  as  dramatic, 
because  the  wavelength  of  the  longitudinal  waves  is  longer  than  that  of  the 
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transverse,  but  further  improvement  is  possible  by  removing  the  discon¬ 
tinuity  in  the  slope^of  u(z)  at  ±t. 

To  be  specific,  Table  II  shows  the  formulas  for  U(k)/U(0)  when 
the  lowest-order  discontinuity  in  u(z)  occurs  in  its  zeroth,  first  and  second 
derivative,  and  also  their  asymptotic  form. 

Table  II 


Three  Fourier  Transform  Pairs 


Since  |k  I  is  of  the  order  of  0.  15  (see  Appendix  I)  and  2 1  might 
be  as  much  as  300  m,  Ik^fl  is  of  the  order  of  i0  or  20.  Thus  equation  (41), 
together  with  the  asymptotic  form  given  for  U(k)/U(0)  in  Table  II,  shows 
that  reductions  of  the  magnetostrictive  ENF  of  tens  of  decibels  are  possible 
by  smoothing  the  sensitivity  profile. 

The  noise  can  also  be  reduced  by  reducing  the  cable  vibration 
level  Sx(<u).  Using  a  mechanically  more  absorptive  cable  or  using  vibration 
isolators  in  the  cable  between  the  antenna  and  the  source  of  the  vibration  are 


ways  of  achieving  this.  If  this  procedure  is  pursued  far  enough,  the  longi- 
tudinal  waves  in  the  antenna  region  will  then  be  generated  locally  by  the 
force  fluctuations  in  the  turbulent  boundary  layer,  as  is  already  the  case 
for  transverse  waves.  Using  the  method  of  the  previous  Sections  one  can 
show  that  the  magnetostrictive  ENF  is  then  given  by 

Seto)-i(^l)  St(0.„)  ±  f  |H{(k,a,)!2dk 


where  S  (k,iu)  is  the  double-spectral  density  of  the  longitudinal  shear  forces 
s 

on  the  cable  and  H^(k,uj)  is  the  spectral  transfer  function  relating  axial,  mo¬ 
tion  to  shear  force  (see  Appendices  D  and  E). 

The  factors  in  the  integrand  here  have  the  same  general  fea¬ 
tures  as  the  corresponding  ones  do  in  the  case  of  transverse  waves  (see 
Fig.  7).  However,  the  wavelength  of  the  longlttltfl/ja)  waves  is  so  much 

i  i2 

longer  than  that  of  the  transverse  waves  that  the  peak  in  the  lH^(k,0))| 

l 

function  occurs  at  a  much  smaller  value  of  k.  The  magnitude  of  |U(k)| 
is  therefore  in  general  bigger  there.  Thus  the  contribution  to  the  integral 

i  1 2 

from  the  peak  in  |H^(k,U)l  is  much  greater  than  that  from  the  peak  in 
2 

|U(k)j  .  The  equation  can  therefore  be  approximated  as 


ro>  B  6  e 


S  (tr)  » 
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S  (0,u>) 
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kfU(kf) 


Ih^u))]2  dk  (43) 


where  k^,  are  the  values  of  k  at  which  the  peaks  of  |H^(k,<«)|  occur.  This 
is  more  conveniently  expressed  in  terms  of  the  mechanical  Q->factor  for 
longitudinal  cable  vibration,  defined  as 
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Then,  since  H^O.w)  =  l/(ou  m^),  where  is  the  effective  longitudinal 
vibrating  mass  per  unit  length  of  cable,  (43)  can  be  rewritten  as 


Sefol  ~  Ef 
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This  formula  is  very  similar  to  (41),  but  now  the  driving  forces  and  cable 
properties  appear  explicitly,  A  general  expression  for  in  closed  form 
is  given  in  Appendix  E  and  a  graph  of  Q^,  for  ^sptcimen  of  the  standard 
0.65  "D  foamed-polyethylene  buoyant  antenna  cable,  is  presented  in  Appen¬ 
dix  I,  together  with  a  similar  graph  of  k^. 

Equation  (45)  shows  that  if  S  (0,u>)  is  known,  then  the  magneto- 
strictive  ENF  of  an  H-field  antenna  can  be  estimated  from  mechanical  and 
electrical  laboratory  experiments  performed  on  the  antenna  [l0].  Since 
the  shear  forces  on  the  cable  depend  only  on  its  diameter  (assuming  it  has 
a  smooth  enough  surface),  then  the  ENF  estimate  obtained  in  this  way 
would  be  a  lower  bound.  Disturbances  propagating  into  the  antenna  region 
(from  sources  oS  vibration  associated  with  the  towing  vessel)  would  gener¬ 
ate  additional  ENF,  but  could  in  principle  be  made  negligible  by  using  a 
long  enough  cable,  using  a  mechanically  more  lossy  cable  or  by  installing  a 
mechanical  isolator  forward  of  the  antenna  region. 

Equations  (41)  and  (45)  suggest  other  means  of  reducing  the 
magnetostrictive  ENF.  The  bias  flux  density  Bq  can  be  minimized  by 
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driving  a  degaussing  current  through  the  antenna  winding  [4]  and  ®  can  be 
reduced  both  by  careful  selection  of  the  core  material  and  by  constructing 
the  core  in  such  a  manner  that  the  core  strain  is  much  less  than-  the  cable 
strain.  Increasing  the  diameter  is  also  a  possibility.  Assuming  that ^ 

S  (0,u>)  has  the  same  behavior  with  radius  as  does  S:(0,oo)  (see  Appendix  D), 
onewcan  deduce  from  (45)  that  the  magnetostrictive  ENF  power  is  in¬ 
versely  proportional  to  radius  squared,  other  things  remaining  constant.^ 

*  »*' 

A  somewhat  more  rapid  reduction  with  radius  than  this  would  be  expected, 
however,  because  the  strength  in  tension  of  the  cable  should  be  proportional 
to  the  radius  of  the  cable.  Thus  the  restoring  force  for  mechanical  vibra¬ 
tion  is  proportional  to  a  ,  but  the  mass  to  be  restored,  per  unit  length*  is 
proportional  to  a  .  Thus  is  proportional  to  (a  /ar  or  a  ,  and  since  for 

a  parabolic  sensitivity  distribution,  k  U(k)/XJ,('0)  is  asymptotically  propor- 
- 1  .  -4 

tional  to  k  ,  then  an  additional  factor  of  a  B  appears  in  the  behavior  with, 
radius  of  the  magnetostrictive  ENF.  A  stiffer,  more  lossy  cable  would 
also  result  in  a  lower  magnetostrictive  ENF. 

The  possibility  of  using  a  strain  relief  core  construction  de¬ 
serves  a  little  farther  attention  in  the  case  in  which  the  core  is  a  helical 
wrap  of  ferromagnetic  tape  By  considering  the  deformation  of  a  helical 
line  in  a  homogeneous  circular  cylinder  a s  rttie  cylinder  stretches,  one 

finds  that  the  strain  6X/£of  the  helical  line  (i.  e.  ,  increase  in  length 
per  unit  length  measured  along  the  helix)  is  given  by 

-j-  -  f  1-  (1  +u)  sin^9  ]  ~  (46) 
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where  u  is  Poisson's  ratio  for  the  material  of  the  cylinder,  0  is  the  helix 
angle  and  6z/z  is  the  axial  strain  of  the  cylinder.  Thus  when  the  helix 
angle  is  zero  (which  means  the  helical  line  degenerates  into  a  straighA^Aine 
parallel -to  the  axis),  the  line  and  cylinder  strain  by  the  same  amount.  If, 
on  the  other  hand,  the  helix  angle  is  ar.c.§vn  (1  +u)”^,  then  the  helical  line 
maintains  a  constant  length  as  the  cylinder  strains.  For  incompressible 
materials,  u  =  1/2,  and  the  corresponding  "zero-strain"  helix  angle  is 


arcsin  \f2/3  or  about  55°. 

Whether  it  is  desirable  to  make  use  of  this  "zero-strain"  helix 

condition  depends  upon  the  relative  levels  of  the  magnetoBtrictive  and  therm- 

.4  ^ 

al  ENF's,.  The  thermal  ENF  power  is  proportional  to  (cos0)  ,  for  a  given 
core  weight  per  unit  length.  This  means  that  taking  0  equal  to  55°  rather 
than  0°  increases  the  thermal  ENF  by  nearly  10  dB. 

In  practice,  the  cere-tape  is  not  an  infinitesimally  thin  Hne  but 
has  a  finite  cross-sectional  area.  This  means  that  the  helix  angle  varies 
across  the  radial  thickness  of  the  core-tape  (since  tan0  =  2iTr/p,  where  r 
is  the  radius  and  p  the  pitch)  so  that  there  will  be  a  varying  linear  strain 


across  the* radial  thickness.  Presumably,  by  arranging  for  0  to  assume  its 
"zero-strain"  value  at  the  radius  of  the  midpoint  of  the  tape  thickness,  the 
net  strain  can  be  minimized.  Other  effects,  such  as  shear  strain  across 
the  width  of  the  core-tape  and  helix  distortion  due  to  the  difference  inme- 
chanical  stiffness  between  the  core  material  and  its  surroundings,  may  also 


need  to  be  taken  into  account. 


The  possibility  of  magnetostrictive  noise  being  generated  by 
transverse  mechanical  waves  should-not  go  without  mention.  However,  if 
the  geomagnetic  field  is  wholly  axial,  the  geometry  has  rotational  symmetry. 
Then*,  the  noise  voltage  generated  by  bending  the  antenna  would  be 
independent  of  the  azimuthal  direction  in  which  the  bend  occurred.  The  voltage 
cannot  therefore  be  linearly  related  to  the  transverse  displacement,  but 
must  depend  on  its  square  or  some  higher  even  power.  It  is  therefore  of 
at  least  second  order,  and  since  the  transverse  displacement  is  very  small 
to  begin  with,  it  should  be  negligible  compared  with  the  longitudinally  gener- 
ated  magnetostrictive  noise.  The  geometry  is  no  longer  rotationaliy 
symmetric  if  a  transverse  component  of  the  geomagnetic  field  exists.  In 
that  case,  however,  the  core  geometry  is  unfavorable  for  the  induction  of 
any  substantial  bias  flux  density  into  the  core,  and  what  does  exist  is  on 
balance,  purely  transverse  (the  axial  components  of  the  induced  flux  density 
are  in  opposite  directions  on  opposite  sides  oC  fhe  core,  if  the  core  is  helically 
wrapped).  Thus  the  only  magnetostrictive  effect  of  bending  is  a  change  in 
the  transverse  flux  density,  which  does  not^c'ouple  with  the  winding.  Another 
rotationaLLy  antisymmetric  geometry  is  one  in  which  the  geomagnetic  field 
is  wholly  axial  but  the  antenna  axis  does  not  coincide  with  the  neutral  axis 
of  the  whole  cable.  Then,  being  offset  from  the  no-strain  axis,  the  core  is 
subject  to  axial  strain  as  a  result  of  puie  bending.  It  is  probably  important, 
therefore,  to  keep  the  axes  of  the  antenna  and  of  the  whole  cable  coincident. 
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VI.  Barkhausen  Noise 

} 

As  the  towing  vessel  changes  its  heading,  the 'smoothly  changing 
axial  component  of  the  geomagnetic  field  is  not  accompanied  by^a/T  equally 
smooth  variation  in  the  resulting  flux  density  it  induces  in  the  antenna  core. 
The  flux  density  changes  in  a  sequence  of  abrupt  jumps  known  as  Barkhausen 
discontinuities  p*  524],  which,  in  turn,  induce  a  wide- spectrum  noise 
voltage  in  the  antenna  winding. 

This  noise  is  the  most  difficult  of  the  known  noise  sources  to 


handle  quantitatively.  This  is  because  the  Barkhausen  noise  generated  by 
different  ferromagnetic  materials  bears  no  obvious  relation  permea¬ 

bility,  and  while  Mumetal,  for  example  [8, has  been  measured  to  have 
very  low  Barkhausen  noise  when  soft,  the  no  se  rises  substantially  when  the 
material  has  undergone  some  work  hardening.  This  work  hardening  would 
occur,  in  particular,  in  the  fabrication  of  a  helical  tape  core. 

The  effect  of  core  area  on  Barkhausen  ENF  is  probably  small, 
although  it  has  not  been  investigated  directly.  This  is  because  the  Bark¬ 
hausen  jumps  tend  to  occur  in  clusters  aB  a  macroscale  effect  so  that 

the  bigger  the  core  cross-section,  the  bigger  the  cluster  involved.  Thus, 
the  noise  increases,  but  the  signal  sensitivity  does  also.  The  clustering 
effect  would  not,  however,  extend  far  along  the  length  of  the  core  and  so 
one  would  expect  the  flux  density  variation  to  be  statistically  independent 
at  widely  separated  points  along  the  antenna.  Thus  Barkhausen  ENF  power 
ic  expected  to  be  proportional  to  the  inverse  first  power  of  antenna  length. 
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A  direct  method  of  preventing  Baridnatisen  noise  is  to  maintain 
the  induced  flux  density  in  the  antenna  core  as  constant  as  possible.  Thus 
an  automatically  controlled  degaussing  current,  which  maintains  antenna 
sensitivity  and  reduces  magnetostr/^/Ve  noise  can  also  reduce  the  Bark- 
hausen  noise.  Measurements  of  the  Barkhausen  noise  caused  by  applied 
field  changes  of  different  rates  enable  one  to  specify  the  performance  re¬ 
quired  from  the  control  loop  [4].  For  noisy  materials,  the  open  loop  gain 
of  the  control  loop  may  turn  out  to  be  unrealistically  large. 

VII.  Some  Experimental  Results 

In  Appendix  D,  a  comparison  is  made  between  the  measured 
noise  of  a  long  E-field  antenna  ancTfne  motion-induced  noise  one  would 
expect  to  measure  if  the  antenna  motion  were  due  to  the  pressure  fluctua¬ 
tions  of  the  turbulent  boundary  layer.  It  was  found  that  the  predicted  and 
measured  noises  follow  the  same  law  with  both  frequency  and  speed 
(namely,  the  noise  power  spectral  densityis  proportional  to  the  fifth  power 
of  towing  speed  and  inverse  fourth  power  of  freauency).  On  the«*o£her  hand, 
there  is  a  constant  difference  between  the  predicted  and  measured  level  of 
some  16  dB.  This  may  be  due  to  the  fact  that  the  pressure  fluctuation  data 
used  in  the  estimate  came  from  the  literature  in  which  only  flat  or  near- 
flat  surfaces  have  been  dealt  with. 

Further  evidence  that  the  transverse  mechanical  waves  on  the 
cable  arise  from  the  turbulence  in  the  boundary  layer  comes  from  measure¬ 
ments  made  with  towed  cables  equipped  with  strain  gages  and  accelerometers, 
It  was  found  that  the  cable  curvature  (measured  with  diametrically  opposed 
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strain  gages)  and  acceleration  at  different  points  along  the  cable  were 
consistent  with  the  hypothesis  that  a  statistically  uniform  driving  mechan¬ 
ism  existed  over  the  whole  cable  length.  The  only  other  possibility,  that 
the  cable  motion  is  excited  near  the  tow-point  or  near  the  free-end  and  then 
propagates  as  transverse  waves  to  the  rest  of  the  cable,  cannot  be  recon¬ 
ciled  with  the  measurements.  This  is  because  the  attenuation  of  transverse 
waves  is  so  great  {see  Appendix^^nat  substantially  different  vibration 

'"V^SSr  ' 

levels  would  be  observed  by  two  sensors  spaced,  say,  a  couple  of  hundred 
feet  apart.  No  such  differences  were  observed. 

For  longitudinal  waves,  the  situation  is  less  clear.  Shear-stress 
fluctuations  in  turbulent  boundary  layers  have  not  apparently  been  treated' 
in  the  open  literature,  so  that  no  predictions  can  be  made  of  the  role  they 
play.  In  addition,  the  attenuation  of  longitudinal  waves  is  small  enough  (see 
Appendix  I)  that  the  effect  of  tow-Tjovnt  disturbances  may  well  be  observable 

.-rr*- 

along  the  whole  cable  length. 

However,  data  obtained  with  strain  gages  show  two  separate 
phenomena  which  suggest  that  local  (i.  e.  ,  boundary-layer)  excitation  is,  in 
fact,  what  predominates.  The  first  phenomenon  is  the  observation  of  a 
marked  standing-wave  pattern  near  the  tow  point.  This  means  that  the 
en&X'gy  is  not  entering  the  cable  at  the  tow  point,  for  if  this  were  the  case, 
the  energy  would  be  flowing  unidirectionally  away  from  the  tow  point  and  no 
standing  waves  would  be  observed.  Thus  the  energy  is  generated  on  the 
cable  somewhere  else  and  then  propagates  along  to  the  tow  point  where  it 
is  reflected  back  again  to  interfere  with  the  incident  energy  and  form 
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standing  waves.  The  second  phenomenon  is  that  the  strain  amplitude  in  the 
cable  is  essentially  uniform  over  its  whole  length.  Thus  the  exciting^source 
must  be  well  distributed  and  statistically  uniform.  (It  is  true  that  the  atten¬ 
uation  per  meter  of  longitudinal  waves  is  small,  but  the  cable  is  so  long  that 
an  easily  measurable  difference  in  level  would  be  present  if  the  energy  were 
generated  at  one  location. ) 

The  strain-gage  data,  together  with  the  theory  developed  in  the 
previous  sections,  allow  one  to  predict  the  noise  voltages  generated  by  the 
different  antennas  tested  experimentally.  Considering  first  motion^induced 

noise,1  since  cable  curvature  is  the  second  derivative  with  respect  to 

•jmi. 


distance  of  cable  displacement,  then  the  curvature  spectral  density 
is  gly„eff  closely  by  S^(uj)/kt^,  where  S^(tu)  is  the  displacement  spectral 
density.  Therefore,  combining  (G6)  with  (E 12 )  one  finds  for  the  motion- 
induced  voltage  spectrum  Sv(<d)  of  a  long  E-field  antenna 


sv(u»  =  !efc  sem  = 
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where  kfc  and  Qt  are  to  be  evaluated  using  the  cable  tension  appropriate  for 
the  position  along  the  cable  at  which  the  curvature  is  measured. 

Similarly,  by  combining  (35)  with  (E12),  and  using  the  results 
of  Appendix  H,  one  finds  for  the  motion-induced  voltage  spectrum  of  a  long 
H-field  antenna  with  a  uniform  profile 
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where  the  suffix  l,  -i  or  g  means  that  the  enclosed  quantity  is  to  be  evalua¬ 
ted  using  the  tension  appropriate  for  the  position  along  the  cable  of  the  two 


ends  of  the  antenna  (£<&  -£)  or  of  the  strain  gages  measuring  the  curvature  (g) . 

In  Fig.  13,  graph  a)  shows  the  cu.vature  spectrum  measured 
with  a  strain-gage  pair  placed  400  feet  from  the  free  end  of  a  cable /being 
towed  at  12  knots.  In  graph  b),  the  result  of  applying  (47)  to  this  curvature 
spectrum  to  derive  the  motion-induced  noise  voltage  of  a  300-meter  long 

* 

n 

E-field  antenna  towed -at  12  knots  is  compared  with  the  "measured'1  E-field 
antenna  noise  voltage.  The.l'^^-Mired 1 1  12  knot  curve  was  actually  ob¬ 
tained  by  lowering  the  14  knot  data  of  Fig.  3  by  3.  35  dB,  since  the  noise 
voltage  follows  a  speed-to-the-fifth-power  law.  Graph  c)  shows  the  result 
of  applying  (48)  to  the  curvature  spectrum  of  graph  a)  to  derive  the  motion- 
induced  noise  voltage  of  an  800-foot  long  H-field  antenna  having  a  uniform 
sensitivity  profile  and  an  effective  length'  of  5.  3  meters  at  45  Hz.  The 
measured  curve  was  obtained  with  the  RCA  Helix  III  antenna.  For  both 
derivations,  the  values  of  k.  and  were  obtained  from  Figs.  8,  9  and  12, 

>  .  i’.M  u  n.,  . 


By  was  assumed  to  be  0.  5  x  io  weber  m  (0.  5  gauss)  and  the  skin-depth 
was  calculated  using  an  ocean  conductivity  of  4  mho/m. 

It  should  be  noted  first  that  the  calibration  of  the  strain  gages 
was  calculated  solely  from  their  nominal  off-  .at  from  the  c«*bie  axis.  In 
fact,  they  are  somewhat  less  sensitive  than  such  a  procedure  would  indi¬ 
cate  because  the  physical  modification  necessary  for  installation  stiffens 
the  cable  locally.  Also,  the  RC A  antenna  cable  was  structurally  different 
from  the  cable  having  the  measured  properties  summarized  in  Appendix  I. 
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Fig.  9.  Mechanical  Q- factor  Qt  vs  frequency  and  tension  for 
transverse  waves  of  the  0.  65  inch  O.  D*  buoyant  cable. 
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Fig.  12.  Cable  tension  as  a  function  of  distance  from  free  end 
and  speed  for  0.  65  inch  O.  D.  buoyant  cable. 
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In  view  of  these  differences,  therefore,  the  agreement  between  derived  and 
measured  noises  in  graphs  b)  and  cl  is.  not  bad.  Until  thermal  noise  takes 
over  at  higher  frequencies,  the  E-field  antenna  noise  level  is  consistent 
with  the  hypothesis  that  it  is  motion-induced,  and  the  same  applies  to  the 
H-field  antenna,  except  that  it  appears  to  b.^Jhe, profile  roughness  contri¬ 
bution  which  takes  over  from  the  long-profile  term  at  higher  frequencies. 

No  sensitivity-profile  data  is  available  for  the  RCA  antenna  and  so  no  esti¬ 
mate  can  be  made  of  the  profile  roughness  contribution. 

The  remaining  two  graphs  in  Fig.  13  concern  longitudinal  waves 
and  their  generation  of  magnetostrictive  noise.  Graph  d)  is  the  longitudinal 
strain  spectrum  measured  with  the  same  two  gages  us,?d  to  obtain  graph  a), 
and  graph  e)  is  the  corresponding  predictionfO|flhe  magnetostrictive  noise 
voltage  spectrum  together  with  the  actual  measured  spectrum  for  the  Lincoln 
Laboratory  50-foot  long  H-field  antenna. 

,  The  prediction  was  made  in  the  following  way.  The  magneto¬ 
strictive  noise  voltage  spectrum  Sv(to)  is  given,  from  (411  as 


S  (to)  =  *  S  (u>) 
v  e  e 


=  T'  ) 


k  U(-k  ) 
S  (to)  —2 - 2. 

x  U(0) 


Jjj  which,  since  the  Lincoln  Laboratory  antennas  have  nominally  parabolic- 
profiles,  U(k)  is  given  by  the  second  line  in  Table  II  in  Section  V.  But  at 
the  same  strongly  biased  condition  (20  mA)  used  for  measuring  the  antenna 
noise,  the  strain  sensitivity  of  a  2 -foot  long,  single-core-bundle  antenna 
element  was  found  to  be  [10]  32  dB//l  v/unit  strain  at  45  Hz.  Since  the 
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effective  length  of  the  antenna  element  is  ie/50,  then  equation  (49)  can  be 


written  for  the  element- as 


»3 2"  =  10  log. 


,  90tt  B  6e  t 
1  o  e 


which  establishes  the  value  of  the  quotient  B  6e  £  /u  in  (49).  (Strain  is 

O  6  cl 

the  length  derivative  of  displacement,  and  |k  |  is  essentially  equal  to  so 

.  ,2  ° 
that  I kQ |  Sx(<u)  can  be' replaced  by  the  strain  spectrum  S^(rn),  and  for  the 

element,  for  which  k^£  <  <  1,  |  U  (-kQ)/U(0)  j is  essentially  equal  to  unity.  ) 

Using  i)  the  value  for  B  6e  i  /u  established  by  (50),  ii)  the 

2  oca  .>*=' 

strain  spectrum  |k  |  S  (u) )  from  graph  d)  and  iii)  the  formula  for  U(k) 
given  in  the  second  line  in  Table  II  in  Section  V,  one  obtains  from  formula 
(49)  the  "derived"  curve  in  graph  e).  To  make  the  null  in  the  "derived" 
curve  coincide  with  the  minimum  of  the  measured  curve,  a  value  for  the 
longitudinal  wave  speed  was  used  which  is  77%  of  that  implied  by  the  data 
given  in  Appendix  I.  This  is  not  unreasonable,  because  at  low  tensions  it 
was  found  that  the  cable  is  less  stiff  in  stretching. 

ThS  remarkably  close  agreement  between  the  "derived"  and 
measured  curves  gives  strong  support  to  the  theory  of  magnetostrictive 
noise  generation. 

No  Barkhausen  noise  has  apparently  been  observed  under  tow‘3 Tig 
conditions,  because  other  noises  have  always  been  dominant.  Barkhausen 
noise  data  measured  in  the  laboratory  is  available,  however,  in  Ref.  [4], 
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VIII.  Conclusion 


The  previous  sections  have  shoim  th'at  an  air-cored  towed 
H-field  antffSna  would  he  inconveniently  large  if  its  thermal  ENF  is  to  be 
low  enough  for  it  to  approach  the  perforsxag^&Qf  the  towed  E-field  antenna. 
In  addition,  an  H-field  antenna  in  the  existing  buoyant  antenna  cable  would 
generate  too  much  motion-induced  ENF  if  it  had  a  uniform  sensitivity  pro¬ 
file.  Thus  a  ferromagnetic  ^ore  and  a  shaped  sensitivity  profile  appear  to 
be  necessary  features.  The  introduction  of  the  core,  however,  adds  two 
more  noise  sources  --  magnetos trictive  and  Barkhausen. 

^The  theory  developed  to  examine  the  various  noise  sources  has 
been  shown  capable  of  predicting  the  noise^j^ages  by  combining  cable  vi¬ 
bration  data  with  laboratory  measurements  on  antenna  samples.  One 
particularly  troublesome  noise  source  has  been  revealed  to  be  motion- 
induced  noise  contributed  by  "profile- roughness". 

Since  the  available  cable  vibration  data  and  theory  indicate  that 
all  cable  motion  is  excited  by  the  stress  fluctuations  in  the  turbulent  bound¬ 
ary  layer,  i£*Jft)uld  appear  that  using  vibration  isolators  or  towing  further 

behind  the  boat  would  be  ineffective  in  reducing  cable  vibration.  Measures 

*■  • 

that  would  be  effective  include  using  a  stiffer,  more  mechanically  lossy 
cable  and  increasing  its  diameter.  Reducing  the  sensitivity  of  the  antenna 
to  the  vibration  is  accomplished  by  tapering  its  profile,  making  its  profile 
as  smooth  as  possible,  using  a  degaussing  current,  using  a  stress- 
relieving  core  construction  and  choosing  a  core  material  low  in  magneto¬ 
striction  and  Barkhausen  noise. 
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APPENDIX  A 

IMML?SION  EFFECT  ON  ANTENNA  IMPEDANCE 

'•■IS8JSK 

(g  j 

An  infinitely  long  solenoid  carrying  a  current  I  and  immersed 
in  a  conducting  medium  induces  in  the  medium  a  magnetic  field  given  by 

H  =  -  ig-  HQ(n(kr)I(m)  (A  - 1) 

at  a  distance  r  from  the  solenoid  axis.  Here,  k  is  the  wave  number  for 

propagation  in  the  medium,  defined  as  yitu  a  pQ  ,  and  I  is  the  axial 

magnetic  current  equivalent  to  the  energized  solenoid,  given  by 
(mi  (el 

I'  =  -  iui^ANI  ,  where  pc  and  A  are  the  permeability  and  cross- 

(VY\\ 

sectional  area  of  the  core.  It  is  assumed  that  I'  is  uniform.  (This 
formula  is  obtained  by  applying  the  principle  of  duality  to  the  well-known 
[2]  result  E  =  -tu  pQH  (kr'  I^6'/ 4  for  the  electric  field  due  to  an  axial 
electric  current.  1 

At  any  radius  r,  this  field  is  the  field  arising  from  all  the  cylin¬ 
drical  shells  of  current  induced  in  the  medium  having  radii  larger  than  r. 
Thus  the  magnetic  field  within  the  solenoid  arising  from  the  external  in¬ 
duced  currents  is  uniform  and  given  by  setting  r  in  (A.  1)  equal  to  r^,  the 
outside  radius  of  whatever  nonconducting  jacket  covers  the  solenoid. 

Hence  the  voltage  AV  induced  per  unit  length  in  the  solenoid  by 
the  external  currents  is  given  by  -itsp^ANH,  and  the  change  in  imped- 
ance  per  unit  length  by  AV/I  .  That  is, 

_  i  2  2  ,  2  ..2  TT  (11,, 

aZ  =  — t  a)  u  ANcrH  'kr  1 
4  rc  o  ,.v  m 
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The  inductance  per  unit  length  hQ  is  easily  shown  to  be  equal  to.N  ^A,  so 
that  the  impedance  change,  normalized  with  respect  to  the  imaginary  com¬ 
ponent  of  the  self  impedance,  is 


=  T-  AH  (1)(kr  )  . 

U)L  4  c  o  m 
o 

Since  the  skin  depth  in  the  ocean  at  ELF  is  many  meters  whereas  r  is 

m 

limited  to  one  or  two  centimeters,' ^i^vSmall-argument  form  of  the  Hankel 
function  can  be  used  to  rewrite  this  as 


AZ 

IB  L 


-  !jl  a.  i.r 

1*0  62  *  L 


in 


VT  j 


-  .577 


in 


+  i7"J 


(A-2) 


At  45  Hz  and  for  the  t/pical  value  of  a  of  4  mho/meter,  the  skin 

-4 

depth  6  is  37.  5  meters.  Thus  if  A  is  one  square  centimeter  or  10  meters, 

p  /p  is  2  x  10^  and  r  is  10”^  meters,  then  AZ/u)  L  is  about  (3.  3  +  i  0.  36) 

co  m  ^ %S$Sfyr-  ° 

x  10  Hence  the  change  in  impedance  is  very  small  compared  with  the 
total  impedance. 

The  thermal  noise  contribution  of  the  water  losses,  expressed 

A 

as  an  equivalent  noise  field  (ENF)  is  given  by  [4  k^  T^  Be  [  AZ}  ]  ®  divided 
by  the  effective  length  per  unit  length  of  the  antenna,  where  k^  is 
Boltzmann's  constant  and  T^  the  temperature  in  degrees  Kelvin.  The 

effective  length  per  unit  length  is,  from  (2),  VTp  AN/(6|a  ),  so  that 

.  r**--  c  ° 


ENF  = 


—  k,  T.  (u  (j.  in  l  — - \  -  . 

’  "  "  °[  V^m) 


-  .  577 


* 


(A- 3) 


is  the  ENF  of  a  1  meter  length  of  the  infinite  solenoid  and  is  proportional 
to  (length) 
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****  For  the  same  assumptions  as  before,  and  for  a  temperature  of 

300°  Kelvin,  this  works  out  to  be  -232  dBE  (dBE  denotes  dB  with  respect 

to  1  volt  per  meter  pei1  V  Hz  ).  For  a  300  meter  length  of  the  solenoid, 

therefore,  the  ENF  contributed  by  the  water  losses  is  -257  dBE.  This  is 

*•*'**' 

lower  than  the  total  ENF  of  the  trailed  300-meter  E-field  antenna  (Fig.  3) 
by  some  50  dB  and  so  is  not  of  concern. 

This  analysis  of  the  effects  of  immersion  has  used  as  a  model 
the  infinitely  long  uniform  solenoid.  It  is  probably  not  accurate  unless 
the  antenna  is  many  skin-depths  long.  However,  since  the  effect  has  been 
shown  to  be  very  small  and  since  the"n2T/?/enna  length  is  likely  to  be  a  few 
skin  depths  long,  one  can  safely  assu^ne  that  the  effect  on  a  real  antenna 
is  not  accurately  known  but  is  still  negligible. 


A^ENDIX  B 

OPTIMUM  AIR-CORE  ANTENNA 

For  an  air-core  antenna,  the  flux  Unking  a  turn  of  radius  r  is 

2  2 
given  by  tt  r  |ioH.,  so  that  the  total  flux  linkage  is  the  quantity  tt  r  (jlqHN  M 

infegrated  over  the  whole  length  of  the  antenna  and  also  over  the  radial 

thickness  of  the  winding.  Here  N is  the  turns  density  per  meter^of  length 

and  M  the  turns  density  per  meter  of  radial  dimension.  Thus,  since 

H  =  Vi^/ E,  the  effective  XengiWtoi  Uie  antenna  is,  in  absolute  magnitude, 

ie  =  17  / NMr2dz  dr  .  (B-l) 


The  resistance  per  turn  is  2-rrrNM/c^,  which  means  that  the  total  resist¬ 


ance  is  given  by 


R  = 


2ir 


w 


/ 


N2M2  *dz  dr  . 


(B-2) 


By  applying  the  calculus  of  variations  to  (B.  1)  and  (B.  2),  one  can  show 
that  the  minirrum  thermal  ENF  power  spectrum,  defined  as  4k^T^R/fe 
is  attained  when  N  is  a  constant  and  M  is  proportional  to  r.  Then  the 
thermal  ENF  power  spectrum  Sg(u))  is  given  by 


o  k,  T,  6 

S  (u;)  =  —  '  -  K- 

e  tr  _  . 

wb4f 


1 


1  -  (a/b> 


-•23 

where  is  Boltzmann's  constant  (1.  38  x  10“  J  j/deg.  ),  T^  is  the  absolute 
temperature  in  degrees  Kelvin,  6  is  *V\e  skin-depth  in  the  ocean  water,  o 
is  the  wire  conductivity,  2f  is  the  total  antenna  length  and  b  and  a  are  the 


(B-3) 
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outer  and  inner  radii  of  the  winding.  (It  is  characteristic  of  all  the  various 
noise-source  ENF's  that  they  are  independent  of  the  turns  density.  Equation 
(B.  3)  exemplifies  this  fact  for  thermal  ENF. ) 

To  obtain  a  lower  bound  for  the  thermal  ENF  of  the  air-core 
antenna,  it  is  convenient  to  assume  that  the  total  cable  weight  is  assigned 
to  thewinding  conductor.  All  the  insulation,  waterproofing  jacket  and  other 
space  within  the  cable  are  assumed  to  be  of  negligible  weight.  It  is  also 
assumed  that  these  materials  are  infinitesimally  thin  where  appropriate. 
Thus  b  is  assumed  to  be  the  outside  cable  radius  and  the  inner  winding 
radius  is  chosen  to  mak2  the  cable  density  less  than  or  equal  to  that  of 
water.  That  is 

(i-  i/p,/. 

0 

where  Pw  is  the  relative  density  of  the  wire  material. 

Thus,  from  (B-3)  and  (B-4),  if  f  =  45  Hz,  a  =  4  mho/m, 

T.  =  300°K,  U  =  300m  and  b  =  0.  325  inches,  then  for  an  aluminum  winding 

n 

(D  =2.7,a  =  3.  53  x  10  mho/m!  the  thermal  ENF  turns  out  to  be  -180.  0 

w  w 

dBE.  The  lighter  (but  impractical)  winding  material  sodium  provides  a 
0.  3  dB  decrease  in  tho  ENF,  while  copper  increases  it  by  2.  2  dB. 


P  >  1 
w 


P  <  1 
w 


(B-4) 
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APPENDIX  C 


SENSITIVITY  DISTRIBUTION  FOR  ANTENNA  WITH. FERROUS  CORE 

If  the  H-field  antenna  has  a  ferromagnetic  core,  the  magnetic 
flux  in  the  core  can  be  calculated  only  by  solving  an  integral  equation. 
Fortunately,  however,  interest  here  is  confined  to  the  total  flux  linkage 
rather  than  the  flux  linkage  per  unit  length.  The  former  is  easily  calcu¬ 
lated  (or  measured)  as  shown  below. 

The  transverse  displacement  of  a  towed  antenna  cable  has  been 
measured  to  be  of  the  order  of  a  micron  (r.  m.  s.  per  root  Herz).  Thus 
the  cable  is  essentially  straight  and  so  the  effect  of  the  flux  in  one  element 
of  the  antenna  core  on  another  an  arbitrary  distance  along  the  cable  from 
it  is  described  by  the  Green's  function  G(z,  zq)  of  an  ideally  straight  core. 

Thus  the  small  transverse  displacement  of  the  antenna  causes 
a  perturbation  of  the  axial  geomagnetic  field  component  of  amount 
(B  /po)5y/dz  and  a  perturbation  of  the  scattered  magnetic  field  of 
/  G(z,so)$  (zQ)dzo,  where  §(z)  is  the  perturbation  of  the  induced  flux 
in  the  core.  The  total  field  perturbation  H  is  the  sum  of  these,  and  at 
the  surface  of  the  core  it  must  be  equal  to  $  times  the  reluctance  per  unit 
length  R  of  the  core.  That  is 

H(z)  =  — ■  +  f  G(z,  z  l  §(z  )  dz  =  R(z)§(z)  . 
u  oz  J  o  o  o 

(C-l) 

A  second  integral  equation  can  be  written  for  the  different 
situation  in  which  the  incident  field  N(z)Iq  is  provided  by  a  current  I  in  the 
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antenna  winding  of  turn  density  N(z).  Then 


Hq(z)  =  N(z)Iq  +  f  G(z,  zQ)  §o(zo)dzQ  =  R(z)  $Q(z)  (C-2) 

where  the  induced  flux  in  the  core  is  denoted  by  $q(z)  and  the  total  field 
by  Hq(z). 

Multiplying  (C-2)  by  $  and  integrating  over  z,  one  obtains 


IQ  y*N(z)$(z)dz  +y*$0(zQ)  /c(zo  ,  z)  $(z)dz"j  dz^  ,  =  Jr(z)  $  q(z)  $(z)  dz, 


(C-3) 


in  which  the  symmetry  of  the  Green's  function  (a  consequence  of  the  recip¬ 
rocity  theorem)  has  been  invoked  to  reverse  the  order  of  the  double  integral. 
When  the  inner  integral  in  (C-3)  is  replaced  by  the  equivalent  expression 
given  by  (C-l),  the  result  is  that  two  of  the  four  integrals  involved  cancel, 
leaving 


dy(z,  t) 

dz 


dz  . 


But  the  left  member  is  just  the  total  change  in  flux  linkage  due  to  the  trans¬ 
verse  deformation  of  the  antenna.  Its  tin  e  derivative  is  therefore  the 
motion-induced  voltage,  v(t).  Hence 


f  $  (z)  >2  . 

v(t>  -  b  / -v-  . 

y  J  u  I  oz  ot 
7  ro  o 


(C-41 


Comparing  (C-4)  with  (11),  one  sees  that  the  equations  are  ident¬ 
ical  provided  A(z)N(z)  of  the  air  core  case  is  now  interpreted  more  generally 

as  the  ratio  $  (z)/p  I  .  But  if  the  length  to  diameter  ratio  of  the  core  is 
o  o  o 

large,  this  ratio  is  simpiy  u?N(z)A(z),  where  is  the  re  rmeability 


t-5 


of  the  core  material.  For  a  composite  core,  pa  is  the  total  flux  within. 
^the  winding  divided  by  the  external  axial  field  strength  and  by  the  permea¬ 
bility  of  free  space 

The  "sensitivity  profile"  $q(z)/|jioIo  can  also  be  directly 
measured  by  passing  an  alternating -current  I  through  the  antenna  winding 
and  measuring  the  flux  in  the  core  from  point  to  point  by  means  of  an  ex¬ 
ternal  coaxial  search  coil. 
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jP*  A-V 


APPENDIX  D 

CABLE  MOTION  DUE  TO  TURBULENT  BOUNDARY  LAYER  STRESSES 

There  are  perpendicular  pressure  fluctuations  and  tangential 
shear  stress  fluctuations  on  the  cable  surface  associated  with  the  turbulent 
boundary  layer  in  the  water  around  the  cable.  There  may  be  additional 
forces  on  the  cable  arising  from  vibrations  of  the  tow  point  and  from  tur¬ 
bulence  induced  by  the  wake  of  the  towing  ship,  but  their  contribution  to 
the  vibration  of  the  cable  in  the  antenna  region  can,  in  principle,  be  made 
arbitrarily  small  by  using  a  long  enough  tow  line.  Thus  a  lower  bound  on 
the  vibration  level  of  the  cable  is  that  generated  by  the  locally-induced 
turbulent  boundary  layer  around  the  cable.  Its  characteristics  are  deter¬ 
mined  onlv  isy  the  towing  speed  and  the  diameter  of  the  cable.  (Presumably, 
the  surface  can  be  made  smooth  enough  that  surface  roughness  is  not  a 
factor.  Franz  [13]  gives  as  a  criterion  of  hydraulic  smoothness  the  con¬ 
dition  that  the  height  of  the  surface  perturbations,  in  inches,  shall  be  not 
greater  than  0.01  divided  by  the  towing  speed  in  knots.  Angle  of  incidence 
seems  also  not  to  be  a  factor  since  the  transverse  drag  component  is  small 
compared  with  the  longitudinal  component  at  the  tow-angles  taken  by  this 
nearly  neutrally-buoyant  cable.  ) 

It  p(9.z,t)  represents  the  pressure  on  the  cable  surface  as  a 
function  of  position  0,  z  and  time  t,  the  transverse  force  per  unit  length, 
assuming  the  shear  stresses  are  relatively  insignificant,  is 

IT 

f(z,t)  =  a  J" p(9,z,t)cos9  d9  • 

-TT 


(>7 


.  .  -*&£- 


The  correlation  function  C^(C,  t)  of  the  transversenj^xte  distribution  is 
given  therefore  by 

Cf(C,T)  =  <f(z,t)f(z+C,  t  +  t )> 

y.  ir 

<p(0,  z,  t)  p(  0 z  +  £,  t  +  t)>  cos  0  cos  0 '  d0  d0 ' 

-ir 

where  the  brackets  <> denote  expected  value.  If  the  transverse  correlation 
distance  is  small  compared  with  the  circumference  of  the  cylinder,  the 
expected  value  in  the  integrand  is  non-negligible  only  when  0  is  close  to  0'. 
In  that  case,  therefore,  C^C.t)  can  be  written  approximately  as 

Cf(c.T)  =  a 1 J  Cp(*,C,T)d*  J  cos20  d0 

-7T  -IT 


(D-l) 


Aleo  implicit  in  this  result  is  the  assumption  that  the  pressure  fluctuations 
are  statistically  uniform  in  azimuth. 

fhe  two-dimensional  spectral  density  S^k.co)  of  the  transverse 
force  distribution  is  simply  the  double  Fourier  transform  of  C^C,  t).  That 
is, 

Sf(k,oo)  =  y'cf{C,T)eiU)T"lkC  dr  dC  , 

where  the  limits  of  integration  are  -  ®  to  os,  or,  from  (D-l) 

ir 

S  (k,u))  =  tra2  f  dTdC  fc  (*,  C,T)eia,T  "ikC  d« 

J  ‘Ctt  P 

-ir 

=  ira2yd<(  /cp(*,C,.)e-kCdC.  (D-2) 
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Here  C,tu^  *s  the  cross-spectral  density  of  the  pressure  fluctuations  at 

two  points  separated  in  azimuth  by  and  axially  by  £,  and  given  by 


Cp(t,C.»)  =  J Cp^.C.Tje^dT  . 

.«*,  It  is,  with  ijt  written  as  ?/a,  the  function  r(to,  5#£j‘1fftrdied  by  Bakewell  [14] 
on  a  truncated  cylinder  of  large  radius. 

There  appear  to  be  no  data  available  on  the  surface  pressure 
fluctuations  on  small  cylinders  and  so  this  derivation  must  continue  on  the 
assumption  that  Bakewell's  data  is  applicable.  The  error  incurred  by  such 
an  assumption  may  be  large,  because  Bakewell  showed  that  his  data  were  in 
very  close  agreement  with  data  obtained  on  flat  surfaces.  The  cable,  however, 
has  a  radius  of  curvature  which  is  certainly  not  large  compared  with  the 
thickness  of  the  boundary  layer,  and  so  cannot  be  regarded  as  flat. 

However,  proceeding  formally,  one  can  express  C  (ijr,C*<D) 

P 

from  Bakewell,  for  -it  &  ^  £  u  ,  as 


Cp(M.<*>)  =  ^(u))A(a)c/Uc)B(u)a^,/Uc)e"lu,^Uc 


(D-3) 


where  Uc  is  the  convection  velocity  at  frequency  u>.  Bakewell  measured 

the  functions  cp ,  A,  B  and  Uc  experimentally.  He  presents  the  results 

graphically,  but  it  would  appear  that  his  A(y)  and  B(|)  functions  are  well 

represented  by  the  expressions  A(y)  =  exp(-(x|/^)  and  B(|)  =  exp(-jf|). 
Substituting  these  expressions  in  (D-3)  and  then  integrating,  in  the  manner 

defined  by  (D-2)  to  obtain  S^(k, u>>,  one  finds 


S f (k, (JJ )  «  a  -4-  <p(ti u) 


0.  01  +  (1  +kU  /to)' 
c 


(D-4) 
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(The  finite  limits  of  integration  in  (D2)  were  extended  to  infinity  to  obtain 
this  result.  This  approximation  is  consistent  with  the  assumption  that  the 
transverse  correlation  distance  is  small  compared  with  the  circumference 
of  the  cylinder.  The  assumption  is,  in  fact,  chi y  marginally  justified  for 
the  correlation  distance  measured  on  flat  surfaces.  It  seems  reasonable 
to  hypothesize,  however,  that  on  this  more  curved  surface  the  correlation 
distance  would  be  smaller. ) 

The  corresponding  spectrum  S^(k, u>)  of  the  cable  displacement 
is  therefore  (see  Appendix  E) 

,2 

S  (k,w)  =  |Ht(k,<u)|  Sf(k,u>)  (D-5) 

where  Ht(k,m)  is  the  transfer  function  for  transverse  displacement,  and  the 
spectrum  S^(tu)  of  the  displacement  of  a  particular  point  is 

V”1'  ■  -h  1 1  Ht(k,t»)|2  Sf(k,u>)dk  .  (D-6) 

By  examining  the  form  taken  by  Ht(k,tu)  (Appendix  E)  and  Sj.(k,u>) 
as  defined  by  (Dr4),  one  concludes  that  for  a  cable  having  properties  similar 
to  those  of  the  existing  buoyant  cable  (Appendix  I),  the  peak  in  Sf(k,u>)  oc- 

2 

curing  at  k  =  -uj/Uc  is  far  enough  away  from  the  region  in  which  iHt(k,iu)| 

has  non-negligible  magnitude,  that  it  is  sufficient  to  assume  S^k.u))  **Sj.(0,u>) 

and  take  it  outside  the  integral  sign.  Since  |  Hfc(k, tu)  |  is  symmetrical,  its 

peaks  at_+k^  pick  up  a  contribution  proportional  to  [S^('k^,^+S^.(k^,  tu)]/2, 

which  is  close  to  being  equal  to  S^O.w)  even  though  S^(-k^,  tu)  >  S^(0,oi)  and 

S.(k  ,  <4  <  S.(0,uij.  Therefore  S  (w)  can  be  written  as 
t  Xj  t  y 
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(D-7) 


»- 


s  (»)  =  Sf(0,(u) 


kfc  ®t 


4,  ,2 

ww  JmJ 


with  kj.,  and  mt  defined  in  Appendix  E. 

The  motion-induced  ENF  of  the  E-field  antenna,  for  example, 
can  be  obtained' immediately  by  combining  (D-4)  with  (G-6).  The  result  is, 
after  some  re-ordering, 


B 


seM 


5ir(l+b ) 


5*. 
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(p(  tu) 
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o2u 


6* 
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u  5 

CO 

TTT 

U)  fa 


(D-8) 

in  which  is  the  towing  speed  and  b  is  the  buoyancy  of  the  cable  (its 
weight  divided  by  the  weight  of  the  displaced  water).  The  length  6#  is  the 
displacement  thickness  of  the  boundary  layer  [15],  and  is  the  thickness 
within  which  most  of  the  change  in  average  water  velocity  between  the  cable 
surface  and  the  water  velocity  at  infinity  has  been  accomplished.  It  will 
probably  be  some  small  par^  of  a  cable  diameter.  The  water  mass  density 
is  denoted  by  p. 

The  factors  in  (D-8)  are  arranged  in  order  of  increasing 
strength  of  variation.  The  first,  involving  and  b,  is  essentially  a 
constant.  The  same  is  true  for  the  second,  except  that  its  value  is  un¬ 
certain.  The  third  and  fourth,  presented  graphically  by  Bakewell  [14], 
are  of  relatively  slow  variation,  while  the  last  is  a  very  strong  function 
of  towing  speed  and  frequency  and  a  fairly  strong  function  of  radius  and 
antenna  length. 
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5  *.  ^ 

The  and  u;  dependencies  appear  to  correspond  fairly  well 
with  the  measured  data .  Figure  3  shows  a  12.  5  dB  increase  in  ENF  as  the 

speed  is  increased  from  8  to  14  knots,  whereas  50  logjQ(14/8)  is  12.2  dB, 

-4 

and  the  slope  of  the  curves  appears  to  follow  the  u>  law.  The  quantitative 
correspondence  is  not  so  good,  however.  Taking  B^  =  0.  5  gauss  =  u 
knots,  b  =  .75,  t  =  150  m,  a  =  0.  325  inches,  f  =  45  Hz  and  5*  =  a/3,  one 
finds  f6*/U  =  0.0172,  so  that  <p(u))/(p2U  ^  iS)  from  Bakewell  [141, 

about  -45  dB  and  Uc/Um  is  0.  87.  This.  leads  to  (tn )  being  -189. 1  dBE 
whereas  the  corresponding  measured  value  (at  14  knots  and  45  Hzl  is,  from 
Fig.  3,  -205  dBE.  Thus  the  estimate  from  turbulence  theory  is  16  dB  too 
large.  This  may  well  be  attributable  directly  to  fhe  fact  that  flat-surface 
pressure  data -is  being  used  to  make  predictions  about  a  strongly  curved 
surface. 

The  derivation  for  longitudinal  vibration  excited  by  axially- 
directed  shear  stress  fluctuations  at  the  cable  surface  is  essentially  the 
same,  with  modified  notation,  except  that  the  cos  0  factor  in  (D-l)  does  not 
appear.  Thus  the  length/ time  correlation  function  C  (£,t)  of  the  axial 
force  per  unit  length  is  given  by 

TT 

Cs(C,T)  =  2TTa2y*Cff(1(,C,T)  d*  ,  (D-9) 

-TT 

where  is  the  correlation  function  of  the  shear  stress  at  points 

separated  in  a*.\77?uth  by  ijt,  length  by  £  and  time  by  t. 

The  double  Fourier  transform  of  C  (£,t)  gives  the  spectral 
density  Sg(k,tu)  of  the  axial  force  cm  i  unit  \ength  <3nd  so  the  spectral  density 
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Sx(k,aj)  of  the  longitudinal  cable  displacement  is  simply 

2 

Sx(k,u>)  =  |  H£(k,cu)| 


(0-10) 


and,  assuming  that  S  (k,u>)  is  essentially  equal  to  S  (0,u>)  over  the  range  of 

S  8  ~ 

,  ,2 

k  for  which  H£(k,u>)|  has  non- negligible  magnitude,  one  finds 


S  (at)  = 
x 


wui4|m£| 


— '  S8(0.») 


(D- 11) 


APPENDIX  E 


CABLE  DYNAMICS 


For  transverse  vibration,  the  cable  is  a  long  beam  with  a 
certain  bending  stiffness  El,  non-uniform  tension  T(z)  and  effective  mass 
(see  Appendix  F)  per  unit  length  excited  into  transverse  vibration  by  a 
transverse  force  distribution  f.  Its  equation  of  motion  is  well  known  [12] 


to  be 


'EI  ,  4  --5T  CT(zl  as 


3  f-n,.,  3  ]  +  -  s 


m 


t  at2 


y(z,  t)  =  f(z,  t). 


(E-l) 

Energy  dissipation  will  limit  the  cable  perturbation  caused 
by  a  single  point  force  at  z  =  zq  to  a  finite  region  containing  the  point  of 
application.  If  the  tension  is  essentially  uniform  over  this  region,  then 
the  cable  perturbation  can  be  calculated  accurately  by  assuming  that  the 
tension  is  truly  uniform  and  has  the  value  T(zq).  Then  (E-l)  can  be  written 


^4  -s2 

EI  ,-V  -  T,V  rr 

3z  dz 


+  m 


t  at2 


y  (z,  t)  =  f  (t)  6(z-z  1 
o  o  o 


(E-2) 


where  f  (t)  is  the  magnitude  of  the  force  and  6(z-z  )  is  the  Dirac  delta 
o  o 

function. 

Taking  the  double  Fourier  transform  of  (E-2),  one  can  rewrite 

it  in  the  spatial  and  temporal  frequency  domain  as  follows: 

Id  ?  ?  \  -ikz 

IEI  k  4  T(z  )k  -  a)  m,  JY  (k,ir)  =  F  (<u)e 

\  °  V  0  ° 
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is  given  by 


Thus  the  spectral  amplitude  YQ(k,o))  of  the  cable  displacement 


Y  (k, (d )  =  H  (k,co)  F  (<u)e  °  (E-3) 

o  to 

where 

Ht(k,«>)  =  - - - - -  (E-4) 

El  k^  +  T(zQ)k2  —  iu^mt 


and  it  is  understood  that  H,(k,<u)  is  also  a  function  of  z  . 

t  o 

Since  Go(z,0))  and  H^k.w)  are  a  transform  pair  in  the  k— •  z 
transformation,  defined  inequations  (24),  then  /  G  (z,u))dz  is  simply  H  (O.iu) 

O  t 

.  2 

which,  from  (E-4)  is  just  -  1/ (<u  mt). 

In  a  manner  entirely  paralleling  that  used  in  the  temporal  fre¬ 
quency  domain  to  study  the  effect  of  linear  filters  on  randomly  time-varying 
input  signals,  one  can  show  that  the  spectral  density  S^(k,tu)  in  the  spatial 
and  temporal  frequency  domain  of  the  transverse  cable  displacement  y(z,t) 
induced  by  the  statistically  uniform  and  stationary  transverse  force  spec¬ 
trum  Sf(k,u))  is  given  by 

S  (k,u>)  =  |Ht(k,u))|2  Sf(k,u>)  (E-6) 

The  resulting  spectral  density  S^,(u)),  in  the  temporal  frequency  domain  only, 
of  the  transverse  displacement  of  a  particular  point  on  the  cable  is  there¬ 


fore  given  by 

Sy(<U)  =  2v  f Sy(k'!i,)dk 

or,  by  virtue  of  (E-6) 


(E-7) 
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(E-8) 


i  /*  2. 

Sy(to)  =  2 v  J  | Hfc(k,  (»‘)*p^(kf  tu )  dk 

In  the  present  application,  Sj.(k,u>)  is  essentially  flat  over  the 
region  in  which. Ht(k,<u)  has  non-negligible  magnitude.  Therefore  (E-8) 
can  be  rewritten  approximately  as 

S  («o)  *  Sf{0,«u)  2|/|Ht(k,u))|2dk  .  (E-9) 

This,  in  turn,  can  be  expressed  in  terms  of  the  mechanical  Q-factor  for 
transverse  vibration,  defined  as 


o  JL  A 

t  '  2  k tJ  |  Ht(0,0J> 


dk  , 


(E- 10) 


where  k^  is  the  (real)  value  of  k  at  which  Ht(k,  u>)  is  a  maximum.  Thus  (E-9) 


becomes 


ktQt  i  ^Z 

S  (u>)  «  |  Ht(0,u>l  j  Sf (0 ,  to > 
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Sf(0,u>)  . 


(E- 12) 


The  peaks  oi  jH^tk,  at  k  =  ±kfc  are  associated  with  poles  of 
the  same  function  lying  near  the  real  axis  at  the  complex  values  of  k  for 
which  the  denominator  on  the  right  in  (E-4)  has  a  zero.  The  zeros  of  the 
denominator  are  the  four  characteristic  wave  numbers  ±kp  given  by 
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(E-13) 
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and,  since  the  poles  near  the  rea';.axis  are  those  at  ±k^,  the  peaks  of 
2 

}Hj.(k,u*)l  occur  at  the  values  of  k  given  accurately  by 

kfc  =  Re{kj}  .  (E- 14) 


By  substituting  (E-4)  into  (E- 10)  and  evaluating  the  integral  by  the  calculus 
of  residues,  one  finds  the  following  explicit  expression  for  Qt 


For  longitudinal  vibratLon,  the  theory  is  simplified  considerably 

because  the  equation  of  motion  is  a  second-order  differential  equation  rather 

than  a  fourth-order  one.  It  is,  explicitly, 

(  *2  2  ) 

AE  ■2y-m,“fir  U(z,t)  --  -s  (as,  t)  (E-16) 

(  dz  *  St*  ) 

where  AE  is  the  stretching  stiffness,  the  effective  vibrating  longitudinal 
mass  per  unit  length,  x(z,t)  the  longitudinal  displacement  and  s  (z,t)  the 
longitudinal  external  force  per  unit  length.  There  is  no  z-dependence  of 
the  coefficients  in  this  equation,  a  further  simplification. 

In  the  two-dimensional  transform  sp^co,  the  equation  becomes 

(AEk2-u>2mf)  X(k,o>)  =  8  (k,«u)  (E- 17) 

so  that  the  transfer  function  H^k,u>)  relating  longitudinal  cable  displacement 
to  externally  apolied  longiiudinal  force  is  given  by 
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Hg  (k,co) 


(E- 18) 


1 

Trn — 2 — 

AE  k  -  u>  m ^ 

Thus  if  the  applied  longitudinal  forces  are  statistically  stationary  and  uni¬ 
form,  characterized  by  a  spectral  density  S  (k,  (t),  the  spectral  density 

S  (k.oj)  of  the  resulting  cable  displacement  is  given  by 

* 

2 

Sx(k,tu)  =  |  {k, to)  j  SB(k,u»)  (E-19) 

and  the  spectral  density  S  (w)  of  the  longitudinal  displacement  at  a  partic¬ 
ular  point  is 

sx(cu)  =  2V/  |H,(k,u,)|2Ss(k,a,)dk  .  (E-20) 

On  the  assumption  that  S  g fk, U) )  is  essentially  flat  over  the 

region  of  the  k-axis  in  which  H.(k,u>)  has  n'.  n-negligible  magnitude,  (E-20) 

•v  * 

can  be  rewritten  as 

sx(®)  *  h  f  lKi(fc*®)|  dkSg(0,U)) 


.  ~c 

"  |  H£(°,W)|  Sg  (0,  tO) 


(E-21) 


kg  Qg 
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Here  is  the  Q-fector  for  longitudinal  cable  vibrations,  given  by 

2 


Qf  =  2  k 
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Hf  (k,co) 


1  IHf(0,<u) 


dk 


(E-22) 
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4  k^Imf  kQl 
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in  which  ±kQ  are  the  poles  of  the  transfer  function  H^(k;u>),  given  by 

(E-23) 

i  l2 

and  ±k^  are  the  locations  of  the  peaks  of  |H^{k,u.')j  ,  given  accurately  by 

kf  =  Re { kQ}  .  (E -24) 
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APPENDIX  F 


EFFECTIVE  VIBRATING  MASS  PER  UNIT  LENGTH 


To  derive  the  form  of  the  effective  vibrating  mass  per 'unit 
length,  the  classical  (laminar)  theory  of  small  vibrations  of  immersed 
objects  will  be  invoked  [11]. 

If  u  =  ;S:u  is  the  transverse  cable  velocity  (at  radian  frequency  u>) 

and  v  is  the  water  velocity,  then  everywhere  but  in  the  (laminar)  boundary 

* v*  ■  >- . 

layer,  V  is  the  gradient  of  a  potential  and  the  normal  component  of  v-r 

must  be  zero  at  the  cable  surface.  Thus  in  cylindrical  coordinates 

2 

lji  (r,  0)  =  -  u  —  cos  0  (F-l) 

«v  2 

v  (r,,  9)  =  7  ^  ~  u  (£  cos  0  +  §  sin  0),  (F-2) 

r 

where  a  is  the  cable  radius. 

The  water  pressure  p  is  given  by  Vp  =  iu>  pv,  so  that 


p(r,0)  =  iu)p^{r,0) 


(F-3) 


Thus  the  perpendicular  force  per  unit  area  on  the  cable  surface  is 
-p(a,9)  in  the  radial  direction,  and  the  tangential  force  per  unit  area  on  the 
cable  surface  is  [v(a,  0)-u]  ‘  $(i  -  i ) o V ou u/ 2  in  the  azimuthal  direction, 
since  fv(a,9)  -  u]  ’  ^  is  the  difference  in  velocity  across  the  boundary  layer. 
Here  P  is  the  water  mass  density  and  u  the  kinematic  viscosity  of  the  water. 

The  total  x-directed  drag  force  on  the  cable  is  therefore 
2tt 

f J  j-p(a»  0)  cos  9  -  [v(a,  0 )  -  u]  ‘  9(1  -  i)o^oiu/2  sin  0  j  ad9  • 

°  (F-4) 


80 


Substituting  (F-l),  (F-2)  and  (F-3)  into  (F-4)  and  carrying  out 
the  integration,  one  finds 

=  -iTa^p[l  +  i(l-i)6  ^/alf-ium)  (F-5) 

^where  6w  =  V^u/cu  and  is  the  skin  depth  for  plane  shear  wave  propagation 
in  the  water. 

If  the  cable  is  being  driven  by  a  uniform  transverse  force  dis¬ 
tribution  f  and  if  mc  is  the  mass  per  unit  length  of  the  cable  alone,  then  the 
total  force  is  f  +  f^  and  by  Newton's  law 

f  +  f d  =  mc(-ia)u), 

•a\uce  -iouu  is  the  acceleration.  From  (F-5),  this  can  be,  written 

f  =  {mc  +Tra2ori  +  i(l  -  i)  6 ^a])  (-^u), 
which  shows  that  the  cable  behaves  as  though  it  had  an  effective  mass  dis¬ 
tribution  mt  (for  transverse  vibration)  given  by 

mt  =  mc  +  ira^  p[l  +  i(l  -  i)  6 ^/a)] 

In  practice,  5 a  <<  1  and  so  this  can  be  expressed  simply  as 

m  =  m  +  ira^  p  (1  +  i  6  /a)  (F-6) 

tew 

Thus  the  cable  behaves  in  transverse  vibration  as  though  it  had 
an  effective  mass  equal  to  the  cable  mass  plus  the  mass  of  the  displaced 
water  and,  in  addition,  has  a  small  imaginary  mass  equal  to  6^/a  times 
the  mass  of  the  displaced  water  lepresenting  the  viscous  energy  dissipation. 

For  longitudinal  vibration,  the  surrounding  water  remains 
essentially  stationary  except  in  the  boundary  layer.  Thus  if  zu  is  the 
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longitudinal  cable  velocity  (at  radian  frequency  10)  then  the  shear  force  per 

unit  area  in  the  z-direction  on  the  cable  surface  is  -u(l  -  i)  pVtuu/2  .  The 

total  z-directed  viscous  force  f,  per  unit  length  of  cable  is  therefore 

d 

£d  =  -Its? 

Following  the  same  derivation  as  before,  one  finds  that  for  longitudinal 
vibration  the  cable  behaves  as  if  it  had  a  mass  per  unit  length  m^  given  by 

m.  =  m  +ira^p(i6  /a).  (F-7) 

few 

It  should  be  noted  that  this  derivation  of  both  mt  and  m^  has 
been  carried  out  using  a  laminar-flow  model.  Intuitively  one  might  expect 
that  the  effect  of  the  turbulence  is  merely  to  superimpose  an  additional 
velocity  field  that  is  moved  bodily  by  the  cable  vibration.  Since  all  dis¬ 
placements  are  small,  the  interaction  between  the  turbulence  and  the  vi¬ 
bration  can  be  assumed  to  be  negligible,  allowing  one  to  combine  linearly 
the  effects  of  each  one  acting  alone.  The  adequacy  of  this  interpretation  is 
not  in  doubt  for  calculating  the  additional  real  part  of  the  mass  in  the  case 
of  transverse  vibration  (see,  for  example,  Wambsganss  and  Boers  [18]).  . 

It  remains  a  hypothesis,  however,  for  the  imaginary  (damping)  part. 


pi(l-i)  6  ^/af-iuju). 
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appendix  g 


MOTION-INDUCED  ENF  IN  THE  E- FIELD  ANTENNA 

The  derivation  of  the  motion-induced  ENF  in  the  E-field  antenna 
follows  very  closely  that  carried  out  for  the  H-field  antenna  in  Section  IV. 

It  is  simpler,  however,  because  the  "profile  roughness"  component  is  not 
present. 

In  place  of  (11) ,  the  fundamental  equation  is 


v(t)  =  B. 


dz  , 


y/«w 

which,  in  the  frequency  domain,  becomes 


(G-l) 


"«5r* 


V(io)  =  -  icu  B_  ^u(z)  X(z,U))  dz 


(G-2) 


which  has  the  same  form  as  (14)  except  that  u(z)  appears  instead  of  u'(z). 
In  these  equations,  u(z)  is  equal  to  1  for  |z|  <  t  and  is  zero  otherwise 
and,  with  the  geomagnetic  Peld  lying  in  the  yz  plane,  only  x-directed  vi¬ 
brations  generate  motion-induced  noise. 

By  following  the  same  derivation  as  before,  but  noting  that  the 
effective  length  is  now  /  u(z)dz  =2  £,  one  finds  the  following  expression  for 
the  equivalent  noise  field  power  spectral  density  SQ (u> ) 


seN 


CD2  B  2 

y 


1  f  1  U(k) 
2tt  ./  IU(0) 


2 

|Ht(k,u>)|  Sf(k,a»)dk 


(G-3) 


in  which 


U(k) 


(z)  e^z  dz 


U 


sin  kl 
kf 


(G-4' 
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It  may  be  shown  by  leaving  u  (  z )  initially  unspecified 

«■  ,  -»• 

that  the  "long  profile"  motion- induced  ENF  of  the  E-field  antenna  is 

2 

minimized  by  the  profile  u(z)  which  minimizes  the  functional  / u  (z)dz/ 

2 

[/u(z)dz]  .  The  calculus  of  variations  then  enables  one  to  verify  that  the 
minimizing  profile  is,  in  fact,  the  one  already  possessed  by  the  E-field 
antenna,  namely,  u(z)  =  1  for  |z|  <  1  ant*  zer0  otherwise. 

The  "shortness"  correction  term  Q (irk^)  in  (G-5 >  is  negli¬ 
gible  for  long  E-field  antennas,  in  which  case  the  motion-induced  ENF  be¬ 


comes  simply 


Se(u» 


■*(- 


u)l  mj 


Sf(0,a>) 


(G-6) 


This  equation  shows  that  the  motion-  nduced  ENF  is  independent  of  any 
mechanical  property  of  the  antenna  cable  except  for  its  diameter  (mfc  is 
essentially  unavailable  for  modification  once  the  cable  diameter  is  fixed, 
since  the  cable  will  invariably  be  close  to  neutrally  buoyant). 

There  is  an  additional  noise  generated  by  the  electrodes  of  the 
E-field  antenna  [6].  It  is  now  believed  that  motion-induced  noise  is  the 
limiting  noi"°  for  the  long  antennas  of  current  interest,  however. 
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APPENDIX  H 


EFFECT  OF  NON-UNIFORM  CA3LL  TENSION 

The  express  :>n  for  the  motion-induced  ENF  of  the  H-field  an¬ 
tenna,  derived  in  Section  IV  and  presented  there  as  equation  (33),  has  three 
terms.  The  first:,  the  "long,  smooth  profile"  contribution  does  not  depend 
on  the  tension  in  the  cable  and  therefore  needs  no  modification.  The  second, 
the  "shortness"  contribution,  arises  from  the  discontinuity  in  the  slope  of 
the  profile  at  each  end  of  the  antenna,  as  shown  by  equation  (26)  through  (28) 
in  Section  IV.  Since  (33)  was  derived  on  the  assumption  of  uniform  tension, 
the  "shortness"  contribution  is  made  up  of  two  equal  amounts,  one  from 

XT'" 

each  end.  The  effect  of  non-uniform  tension,  therefore,  is  to  make  the 
contributions  from  each  end  of  unequal  magnitude.  Thus  3Q^irk^jf)  must 
be  replaced  by  3  {  [Qt(f)/kt(f)]  +  [Qt(-£)/kfc(-J?)3  }/(2tt  £),  since  both  and 
kt  depend  upon  the  tension. 

Finally,  considering  the  "profile  roughness"  contribution,  one 

notes  that  since  this  is  generated  over  the  whole  length  of  the  antenna,  the 

tension-dependent  terms  must  be  averaged  appropriately  over  the  antenna 

length.  This  is  accomplished  by  replacing  the  tension-dependent  factor 

k  ^  Q.S  (k.)  by  the  average  /  u^(z  )k,^Q,  S  (k.)dz  / !  \x"  (z  )dz  .  The 
t  t  g  i  O  t  t  c  c  o  oo 

rationale  for  this  is  that  the  "profile  roughness"  contribution  to  the  mean 

2 

square  noise  voltage  is  proportional  to  the  factors  (i)  k.  S  (k, )  (because  it 

t  g  t 

is  the  slope  of  the  profile  roughness  which  is  the  important  quantity), 

tii)  u^(z)  (because,  from  (29),e(z'  is  the  fractional  perturbation  of  the  profile), 
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and  (iii)  k.  Q  {because,  from  (E-12),  the  local  mean  square  displacement 
1  t 

Sy(a?)  is  proportional  to  the  tension-dependent  quantity  kfcQt  as  well  as  to 
other  tension-independent  ones).  In  addition,  the  vibration  at  widely  sepa¬ 
rated  points  is  uncorrelated.  Therefore,  the  total  profile  roughness '’con¬ 
tribution  to  the  mean  square  motion-induced  noise  is  proportional  to  the 
integral  over  the  antenna  length  of  the  product  of  the  factors  i,  ii  and  iii 


listed  above. 

The  "shortness"  correction  term  Qj/.(irk for  the  E-field 
antenna,  given  in  (G-5),  also  reauires  modification  to  include  the  effect 
of  non-uniform  tension.  Since  the  term  arises  in  the  same  way  as  that 


for  the  H-field  antenna,  namely  from  the  profile  discontinuity  at  each  end, 
then,  the  modification  is  effected  in  the  same  way.  Thus  the  term  becomes 
[Qt(f)/kt(f)  +  Qt(-f)/kt(-^)]  /  (2-irf).  The  "long,  smooth  profile"  term  (the 


"1"  in  (G-5)),  requires  no  modification  since  it  is  tension-independent. 


APPENDIX^ 

MECHANICAL  PROPERTIES  OF  0.65  INCH  O.  D.  BUOYANT  CABLE 


To  provide  a  quantitative  context  for  the  theory  developed  1 
the  main  text,  the  basic  mechanical  characteristics  of  the  existing  0.  65  O.  D. 
buoyant  antenna  cable  are  listed  in  this  Appendix. 

The  particular  cable  specimen  to  which  the  following  data  apply 
is  that  specified  in  Ref.  [16],  It  consists  of  a  core  of  four  24  AWG  insulated 
copper  conductors  in  a  0.  165  inch  O.  D.  poi %jSj$&y\zne  jacket  surrounded  by 
a  layer  of  18  fiberglass  strands  to  provide  the  necessary  tensile  strength 
and  jacketed  overall  by  a  0.65  inch  O.  D.  foamed  polyethylene  outer  layer 
to  provide  the  necessary  buoyancy. 

Laboratory  measurements  at  room  temperature  on  a  piece  of 
this  cable  indicate  that  the  complex  stiffness  of  the  cable  in  bending  (El) 
and  stretching  (EA)  are  essentially  constant  over  the  frequency  band  of 
interest  here  (20  to  200  Hz)  and  have  the  values 

El  =  1.2(1  —  i  0.  12)  Newton  (I -  1 ) 

EA  =  0.63  X  106  (1  -i  0.  025)  Newton. 


(If  the  tension  is  less  than  about  500  Newtons,  the  stretching  stiffness  is 

less  than  this  value,  presumably  due  to  slack  strength  members.  )  The 

cable  radius  is  0.00825  m  and  its  average  relative  density  is  0.75,  and 

3  3 

the  mass  density  and  viscosity  of  ocean  water  can  be  taken  to  be  10  kg/m 

.  in-6  2  -1 
and  10  r .  s 


Using  these  values  to  calculate  the  wave  number,  mechanical 
Q-factor  and  attenuation  for  both  transverse  and  longitudinal  waves,  one 
obtains  the  results  presented  in  Figs.  8  through  11.  The  wave  numbers 
and  O  -factors  are  as  defined  by  (E-14),  24),  (E-15),  (E-22)  in  con- 

junction  with  (E-13),  (E-23),  (F-6)  and  (F-7).  The  attenuations  are  defined 
by 


=  20  ^gjolexpClmfkj}]) 

=  20  log10(exp[lm{ko}])  . 


(1-2) 


Since  the  cable  is  uniform  in  its  properties  and  has  (presumably) 
r>'i  drogue  at  the  far  end,  it  has  a  uniform  drag  per  unit  length  and  uniform 
buoyancy  per  unit  length.  If  the  submarfiie'maintains  a  constant  depth,  di¬ 
rection  and  speed  therefore,  the  cable  takes  the  form  of  a  straight 
line  and  the  tension  increases  linearly  with  distance  from  the  free  end.  From 
Ref,  [17*]  one  can  show  that  the  angle  6  (in  radians)  from  the  horizontal 
is  given  approximately  by 


9  = 


jiunm 

i  1. 1 


(1-3) 


where  b  is  the  relative  density  of  the  cable  and  aits  radius,  g  is  the  acceler¬ 
ation  due  to  gravity  and  U  is  the  towing  speed.  This  can  be  expressed  as 


fio  27 

0  «  -£j — 


(1-4) 


knot 


where  9°  is  the  angle  A  expressed  in  degrees  and  U  j.  *s  the  speed  Uw 
expressed  in  knots. 

Again,  using  Ref.  T17],  one  can  show  that  the  tension  T  is  given 


by 
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